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ABSTRACT 

We have studied the optical spectra of a sample of 31 O- and early B-type stars in the Small Magellanic Cloud, 21 of which are associated 
with the young massive cluster NGC 346. Stellar parameters are determined using an automated fitting method (Mokiem et al. 2005), which 
combines the stellar atmosphere code fastwind (Puis et al. 2005) with the genetic algorithm based optimisation routine pikaia (Charbonneau 
1995). Comparison with predictions of stellar evolution that account for stellar rotation does not result in a unique age, though most stars are 
best represented by an age of 1-3 Myr. The automated method allows for a detailed determination of the projected rotational velocities. The 
present day v T sin i distribution of the 21 dwarf stars in our sample is consistent with an underlying rotational velocity (v T ) distribution that can 
be characterised by a mean velocity of about 160-190 km s _I and an effective half width of 100- 150 km s -1 , The v t distribution must include a 
small percentage of slowly rotating stars. If predictions of the time evolution of the equatorial velocity for massive stars within the environment 
of the SMC are correct (Maeder & Meynet 2001), the young age of the cluster implies that this underlying distribution is representative 
for the initial rotational velocity distribution. The location in the Hertzsprung-Russell diagram of the stars showing helium enrichment is in 
qualitative agreement with evolutionary tracks accounting for rotation, but not for those ignoring v t . The mass loss rates of the SMC objects 
having luminosities of logL^/LQ > 5.4 are in excellent agreement with predictions by Vink, de Koter & Lamers (2001). However, for lower 
luminosity stars the winds are too weak to determine M accurately from the optical spectrum. Two of three spectroscopically classified Vz stars 
from our sample are located close to the theoretical zero age main sequence, as expected. Three additional objects of lower luminosity, which 
are not given this classification, are also found to lie near the ZAMS. We argue that this is related to a temperature effect inhibiting relatively 
cool stars from displaying the spectral features characteristic for the Vz luminosity class. 
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1. Introduction 

Mass loss and rotation play a crucial role in the evolution 
of the most massive stars. These two processes are linked. 
To illustrate this, fast rotation may lead to an enhanced mass 
loss, as is, for instance, o bserved in the case of Be stars 
jLamers & Pauldrachlll99ll) and ^Carinae dSmith et al.ll2003l 
Ivan Boekel et al.ll20 03). Conversely, the stellar winds of these 
stars will cause loss of angular momentum, leading to spin 
down. In a Galactic environment the spin down of massive 
stars is expected to be relatively rapid, essentially "erasing" 



their initial rotatio nal velocity p roperties within only a few 
million years ( Mevnet & Maeder 2000). For an environment 
characteristic of the Small Magellanic Cloud (SMC) it is 
predicted that the initial conditions of rotational velocity re- 
main preserved during t he main sequence life of O-type stars 
dMaeder & Mevne3l200ll) . as the strength of st ellar winds di- 
minish es with d ecreasin g metal cont ent (e.g. lAbbottlll982l 
iKudritzki et alll987tlvink et al.l200ll) . 

The study of a young SMC cluster containing a substantial 
population of O-type main sequence stars would allow the role 
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of metallicity on mass loss and angular momentum loss during 
the early evolution of massive stars to be better constrained. 
This provides a better insight into the current and perhaps even 
initial rotational conditions from which this evolution unfolds. 
Among others, this would provide important constraints on the 
physics of the formation of massive stars, for instance on the 
role of magn etic fields and t he lifetime of accretion disks (see 
e.g. lPorter & Riviniusl2003l) . 

The distribution of projected rotational velocities v r sin i of 
Galact i c early-type s tars has been studied bvlConti & Ebbets 
Jl977h . |Pennvl ( ll996l) and lHowarfh etaD i ll 997b . The first extra- 
galactic studies of rotational velocities have been presented by 
iKellerl J20()4) for a sample of Large Magellanic Cloud (LMC) 
B-ty pe stars in the vi cinity of the main sequence turnoff and 
bv IPennv et alJ J200 4) for predominately giant and supergiant 
SMC and LMC O-type stars. So far, no systematic study for 
early main sequence O-type stars in a low metallicity environ- 
ment has been attempted. 

The young cluster NGC 346 in the Small Magellanic Cloud 
(SMC) contains a substantial population of O-type main se- 
quence stars. I n the context of t he VLT-FLAMES Survey of 
Massive Stars JEvans et alJ we have used the Fibre 

Large Array Multi-Element Spectrograph at the ESO Very 
Large Telescope to obtain spectra of several dozen dwarf O- 
type star s in this cluster. Ou r genetic algorithm based fitting 
method (Mok iem et aHl2005l) is used to perform a homoge- 
neous spectroscopic analysis of this set of stars, including the 
projected rotational velocity as one of the fitting parameters. 
The i> r sin; information obtained from this analysis may be 
compared to models of the present and possibly initial (see 
above) v r sin i distribution. This is a main aim of this paper. 

The strength of radiatively driven winds is predicted to be 
reduced for decreasing metal content. In the last two decades 
this prediction has been quantified by different groups who 
found a met allicity d ependen ce of the mass loss rate of M(Z) oc 
Z - 5 '" - 7 (e.g. lKudritzki et alJll987t IPuls et alJEoOot IVink et all 
l200ll) . Qualitatively this metalli city dependence has been con- 
firmed by several authors ("e.g. IPuls et alJll996t iBouret et"all 
120031 lEvans et al]l2004al iMassev et alJl2005h . However, until 
now a quantitative comparison of the theoretically predicted 
and the empirically determined Z dependence is still lacking. 
The analysis of our SMC targets will provide insight in the 
wind characteristics of objects with a metal content approxi- 
mately five times lower than in Galactic objects. Consequently, 
given the large number of objects that we have analysed in a 
homogeneous way, the current study will for the first time be 
able to provide such a quantitative comparison. 

While comparing mass loss rate predictions with observed 
values it is important to realise that recent studies of the wind 
strengths of not too lumino us (logL^/Lp < 5.3) Galactic 
and SMC dw arf O-type stars JBouret et all2003l [Martins et alJ 
120041 l2005bl) seem to indicate that the mass loss rates are sig- 
nificantly (up to two orders of ma gnitude) lower than expected 
from theory llVink. de Koter. & Lamers l200ll) . If this is indeed 
the case it would imply that our chances of actually observ- 
ing a v, pattern that is representative of the initial distribution 
will increase, as less angular momentum will be carried away 
in the stellar wind. However, from the viewpoint of our under- 



standing of the physics of relatively weak stellar winds (~10 
- lO _6 M0yr _1 ) the situation is obviously less desirable. 

The mas s-loss rates reported by IBouret et alJ and 
iMartins et alJ are based on an analysis of unsaturated ul- 
traviolet resonance lines. The M sensitivity that can be 
obtained using such lines is ~1(T 9 - 10~ 8 Moyr~'. This is, 
in principle, much better than is possible using Ha as the 
mass loss diagnostic. Using our automated fitting method and 
spectra having a signal-to-noise of 50-200 (typical for those 
secured within the context of our VLT-FLAMES programme) 
we can push the Ho- method as low as ~2 x 10~ 7 MQyr~ l . 
This provides some overlap with the weak wind regime so 
far reserved for the UV line method. It thus allows for a first 
investigation of the question whether the so-called "weak 
wind problem" points to missing physics in the theoretical 
predictions of mass loss or that there may be problems with 
the mass loss diagnostic, i.e. that we do not fully understand 
the formation of UV resonance lines. 

The structure of this paper is as follows: in Sect. |2]we de- 
scribe the data set that has been analysed using our genetic al- 
gorithm based fitting method, which is discussed in Sect.|5] The 
stellar properties of our sample determined in this analysis are 
presented in Sect.|4] In sections |5]|6] and we investigate and 
discuss, respectively, the underlying rotational velocity distri- 
bution of our sample, the mass loss rates we have determined 
in the context of weak winds and the evolutionary status of the 
cluster NGC 346. Sect.|8]summarises and lists the conclusions 
of our study. Finally, in the appendix the fits and comments on 
the individual objects are provided. 

2. Data description 

The sample considered here is mainly drawn from the tar- 
gets observed in the SMC as part of the VLT-FLAMES survey 
of massive stars ("see lEvans et al . 2005). The survey observed 
two fields in the SMC, centered on the clusters NGC 346 and 
NGC 330. Here we analyse each of the O-type stars and three 
luminous B-type stars observed in the NGC 346 field, and two 
O-type stars from the older NGC 330 cluster. 

To improve our sampling of luminosity and tempera- 
ture within the O-type domain the FLAMES targets were 
supplemented by eight field star s from the catalogue of 
lAzzooardi & Vigneaul i 19751 Il982l hereafter AzV, note that 
in other studies these are also often denoted by AV). These 
stars were observed using the Ultraviolet and Visual Echelle 
Spectrograph (UVES) at the VLT as part of the programs 67. D- 
0238, 70.D-0164 and 074.D-0109 (P.I. Crowther). 

Table prists some of the basic observational properties of 
the programme stars analysed together with common aliases 
used in other studies. For a full descrip tion of the observa tional 
properties of the FLAMES targets see lEvans et alJd2006h . Note 
that five of these objects and six of the field stars were recently 
also analysed using line blanketed stellar atmosphere models 
dCrowfher et alJl2002t IBouret et alJEooi lEvans et al.ll2004at 
Mas sev et al.l l2004: He ap et alJl2006l) . A comparison with our 
findings is presented in the appendix. Each of the FLAMES 
targets were observed with the Giraffe spectrograph at least six 
times, at each of six wavelength settings. The effective resolv- 
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Tabl e 1. Basic paramete rs. Identificati ons: FLAMES targets ("NGC330" and "NGC346 "1 from lEvans et alJ d2006h . "MPG" 
from lMassev et alJ lfl989). "AzV " fromlAzzopardi & Vigneai] ill 9751 fl982l) and "Sk" fromlSanduleak Jl968lK Photome tric data 
for the FLAM ES targets is fro m lMomanv et al.l d200lh : for objects with a pri mary AzV identific ation these are from|Massev 
J2002I) and lMassevetail 12004). Spectral types of the FLAMES targets are from Evans et al. (2006). For convenience previously 
published spect ral types ar e also g iven. Sources of these classificati ons are W77 JWalborrJll977l): NMC JNiemela et alJll986l) : 
G87 JGarmanv et aljll987l) : MPG dMassev et alJll989t) : WO O JWalborn et alJ|200(A : EH 04 llEvans et alJl2004bh . Spectral types 
for the field O-type stars were t aken from the catalogue s of Garmanv^^J (ll^8jb and Walborn et al. (2002). Win d velocities 
given without brackets are fromlGar manv & Fitznatrickl Jl988l) :lHaser et alJ Jl998h : IPrinia & CrowfheJ dl998h : ICrowfher et alJ 
j2002l) : lBouret et alJJ2003l) : lEvans et al.N2004ald) : lMassev et al J ( 120041) . Values between brackets are calculated from the escape 
velocity at the stellar surface. 
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ing powe r of the observatio ns is R 20 000, full details are 
given by Eva ns et all (|2£)06). The multiple exposures, often at 
different epochs, allowed for the detection of va riable r adial 
veloc ities, with a number of binaries detected (see IE vans et alJ 
120061) . Three of the stars analysed here are detected as binaries 
(NGC 346-001 & NGC 346-025) or are found to have variable 
velocities suggestive of binarity (NGC 346-007), as noted in 
Table [0 Note that none of these appear to have massive com- 
panions, which could have a significant effect on the derived 
parameters if they were present. 

A general description of the reduction of the FLAMES data 
was given bv lEvans et al.l d2005l) . The most pertinent part of 
the reductions for the current study is that of sky subtraction. 
A master sky-spectrum was created from combining the sky 
fibres (typically 15), individually scaled by their relative fi- 
bre throughput. In general the sky background is low for our 
(relatively bright) targets, and this approach successfully re- 
moves the background contribution. However, in regions such 
as NGC 346 which has strong nebular emission, accurate sub- 
traction of nebular features in multi-fibre data is notoriously 
difficult. This does not hamper our analysis, except in the core 
of the Ha Balmer line. For most stars the core nebular emis- 
sion is well-resolved and we simply only consider the wings of 
the profile in the automated line fits. Tests assessing the impact 
of possible residual nebular contamination in the line wings or 
over-subtraction of sky components are discussed in Sect. 13.31 

For each wavelength range the individual sky-subtracted 
spectra were co-added, then normalised using a cubic-spline fit 
to the continuum. These normalised data were finally merged to 
obtain a final spectrum covering 3850^1-750 and 6300-6700 A. 
The combined spectra have a typical signal-to-noise ratio of 
50-200, depending on the magnitude of the target. 

Spectral types for the field O- t ype st ars w ere taken from 
the c atalogues of iGarmanv et alJ dl987l) and IWalborn et alJ 
il2002l) . The FLAMES spectra were classified by visual in- 
spection, using published sta ndards, in particular the atlas of 
IWalborn'& Fitzpat rickl dl990l) . with cons ideration to the lowe r 
metallicity environment of the SMC (e.g. Walb orn et ali20 00). 
The classifications quoted in Table ^ are from lEvans et alJ 
(2006). 

The majority of the field stars were observed with the VLT- 
UVES spectrograph during a visitor run on 27-28 September 
2001 under program 67.D-0238. UVES is a two arm cross- 
dispersed echelle spectrograph, the red arm of which contains 
a mosaic of two detectors. A standard blue setting, with cen- 
tral wavelength 437 nm provided continuous coverage between 
3730-5000 A, recorded on a single 2x4K EEV CCD. A non- 
standard red setting with central wavelength 830 nm included 
an identical EEV CCD covered 6370-8313 A, whilst a 2x4K 
MIT-LL CCD covered 8290-10250 A. A 1" slit was used in 
variable seeing conditions to provide a spectral resolution of 
0.09 A at He. Exposure times ranged from 1200 to 3000 sec- 
onds. 

Subsequently, AzV 388 was observed with UVES on 5 
Dec. 2002 under service program 70.D-0164, using a standard 
setup with central wavelength 390 nm (blue, 3300^1-500 A) and 
564 nm (red, 4620-5600 and 5675-6650 A), plus a second red 
setup with central wavelength 520 nm (4170-5160 and 5230- 



6210 A). The exposure time was 2700 s for each setup, with 
a 1.2" slit. Finally, AzV 95 was observed with UVES under 
service program 074.D-0109 on Nov. 27 2004, using standard 
setups with central wavelengths 390/564 nm and 437/860 nm, 
with individual 2200 s exposures, again with a 1.2" slit. In 
all cases, the two-dimensional CCD frames were transformed 
to extracted one-dimensional spectra using the UVES pipeline 
software. The S/N in individual blue spectra ranged from ~50 
to -120. 

The quoted photometry for the FLAMES targets in the 
NGC 346 and NGC 330 fields is taken from the ESQ Imaging 
Survey (EIS) pre-FLAMES release bv iMomanv et alJ fcOOll)" 
Thes e data co mpare well to the values obtained previously by 
iMassev et all dl989l) and lAzzooardi & Vigneaul dl975l Il982). 
with an average difference of 0.04 magn for both V and (B - V) 
and incidental (two cases) maximum differences of approx- 
imately 0.1 magn. Photometric data fo r the field stars was 
taken from the UB VR CCD survey of IMassev! (|2002) and 
IMassev etalld2004 . 

The interstellar extinction (Ay) l isted in T able [T] was cal- 
culated using intrinsic colours from [johnsonl 1 19661 and ref- 
erences therein) and by assuming a ratio of total to selec- 
tive extinction of Ry = 3.1. Using these Ay values, extinc- 
tion corrected visual magnitudes (Vq) were calculated from 
the observed V-band magnitudes. Finally, the absolute vi- 
sual magnitude My, was calculated tak i ng a distance mod- 
ulus of 18.9 magn dWesterlundl Il997l lHarries et all 120031: 
iHilditch et alJ2005h . 

3. Analysis method 

To analyse the large number of spectra we use an automated 
fitting method. This me thod was developed by Mo kiem et alJ 
(2005, hereafter Paper I) for the fitting of the profiles of hydro- 
gen and helium lines and provides a means for an unbiased and 
homogeneous analysis of large samples of early-type stars. For 
a detailed description we refer to this paper. Here we will only 
give a brief description of the method and modifications of its 
implementation. 

The automated fitting method consists of two main com- 
ponents. Th e first componen t is the stellar atmosphere code 
fastwind of IPuls et alJ d2005l) . This fast performance code in- 
corporates non-LTE and line blanketing following the concept 
of a unified model atmosphere to synthesise hydrogen and he- 
lium line profiles. To optimise the input parameters for fast- 
wind, i.e. to determine the fit parameters for a certain spec- 
trum, the genetic a lgorithm based optimisation routine pikaia 
(Charb onneaull 19951) is used. This second component is capa- 
ble of global optimisation and in combination with fastwind 
provides, as was shown in lPaper \ a robust method for the fit- 
ting of the hydrogen and helium spectra for a broad range of O- 
and early B-type stars. 

In short the automated method determines the best fit by 
calculating consecutive generations of fastwind models. In ev- 
ery generation the models which fit the observed spectrum the 
best are selected and their parameters are crossbred and mu- 
tated to create new sets of parameters. Using these sets a new 
generation of models is calculated. The procedure is repeated 
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until the fit quality of the best-fitting model is maximised. This 
fit quality is defined as the inverted sum of the reduced chi 
squared, values of the observed hydrogen and helium line 
profiles and the synthetic line profiles. In this sum each spec- 
tral line has a weight assigned. These weights are used to ex- 
press the accuracy with which the model atmosphere code is 
believed to be able to reproduce certain lines. For instance, the 
He i line at 447 1 A is given a low weight for late type super- 
giants because of complicat ions due to the so-called "gener- 
alised dilution effect" (Voel s et alJll989l) from which this line 
suffers. A lower weight is assigned to the neutral helium sin- 
glet lines fo r early and mid-type b ecause the codes fastwind 
and cmfgen dHillier & Millerlll998|) show a discrep anc y in the 
predictions of these diagnostics JPuls et alJl2005l) . In IPaper 3 
the full weighting scheme is described and discussed. 

3.1. Fit parameters 

For the fitting of the spectra we allow for seven free parameters. 
These parameters are the effective temperature T e ff, the surface 
gravity g, the helium over hydrogen number density Yu e , the 
microturbulent velocity u tur b, the projected rotational velocity 
v r sin i, the mass loss rate M and the exponent of the beta-type 
velocity law describing the s upersoni c wind regime of the stel- 
lar atmosphere. In contrast to lPaper j v T sin i is now treated as a 
free parameter. This could be done because the FLAMES spec- 
tra are of sufficiently high resolution to allow for a self consis- 
tent determination of this parameter (see Sect. l3.2t . Note that 
this parameter can also be determined relatively accurately us- 
ing alternative methods (see Sect. l4.7l i. The reason that we still 
incorporate it as a free parameter is that it allows for a mean- 
ingful estimate of the error in this parameter. More importantly, 
it allows for the propagation of this uncertainty in the error es- 
timates of the other fit parameters. 

An important parameter which cannot be determined from 
the optical spectra of O stars is the terminal velocity of the 
wind Doo. Therefore, if a value obtained from the analysis of ul- 
traviolet (UV) wind lines is available for a certain object, we 
keep v x fixed at that value. If no ««, determination is avail- 
able, the scaling relation of y M with the escape velocity at the 
stellar surface (v esc ) is used throughout the fitting process. For 
early-type s tars this scaling im plies that the ratio Uco/uesc = 2.6 
is adopted jLamers et alJll995ft . In some cases this produced 
relatively large wind velocities (see Tab.Q. However, these in- 
dividual cases do not have a significant impact on our results, 
as for most of these objects we could only derive upper limits 
for the mass loss rate (also see Sect. l4.6t . 

As the current implementation of the fitting method only 
analyses the hydrogen and helium lines, no explicit abundance 
values other than the ratio of these two elements can be deter- 
mined. Therefore, we adopted fixed values for the atmospheric 
abundances of the background m etals. For these value s we use 
the Solar abundances from Grevesse & Sauvall 11998. and ref- 
erences therein) scaled proportionally with respect to mass ra- 
tios by the same factor for all elements heavier than helium. 
Iron, due to its strong line blanketing effect on the stellar at- 
mosphere and emergent spectrum, can be considered to be the 



most important metal element. Large differences between the 
abundances of other metals, such as nitr ogen, have been re- 
ported for the Galaxy and SMC (e.g. lTrundle et alJl2004 . al- 
though their effect upon the spectroscopic analysis of the hy- 
drogen and helium lines are negligible. Consequently, we set 
the metallicity scaling factor equal to the iron abundance ratio 
of the SMC. From the analysis of early type SMC star s this 
ratio is found to be 1/5 times Solar (Rollest on et all2003h . 

3.2. Formal tests 

As argued in Paper I so-called formal tests, i.e. fitting of syn- 
thetic data, are an integral part of the automated fitting method. 
First of all such tests are necessary to assess whether the data 
quality is sufficient to secure a successful determination of the 
global optimum in parameter space, i.e. whether it will recover 
the global best fit. Secondly, they are necessary to estimate the 
minimum number of generations that have to be calculated in 
order to find this best fit. Consequently, by fitting synthetic data 
with a similar quality as our observed spectra, we can establish 
the minimum number of generations that have to be calculated 
to safeguard that the best possible fit will be obtained when 
fitting the real spectra. 

Similar to IPaper J three data sets were created based on 
fastwind models with parameters representing different types 
of early-type stars. For these three data sets, which are denoted 
by A, B and C, the input parameters of the fastwind models 
are listed in Tab. |2] Set A and set B correspond to bright hot 
supergiants. The parameters in set C represent a cooler dwarf 
O-type star. For t he mass loss rates we adopted values based on 
the prediction of IVink et al.l d200ll) assuming Zsmc =0.2 Zq. 

Using the line profiles from the fastwind models the syn- 
thetic data was created by first convolving the profiles with a ro- 
tational broadening profile. To the broadened profiles Gaussian 
distributed noise corresponding to a signal-to-noise ratio of 50, 
was added. This value approximately corresponds to the low- 
est S/N in our sample. Finally, as nebular emission requires us 
to ignore the cores of hydrogen and neutral helium lines in the 
fitting of our target stars, we also removed these cores from our 
test data set. In case of the He i lines 2 A from the central core 
was cut out. From the hydrogen lines, with exception of Ha, 
3 A was removed. Of all observed line profiles Ha exhibits the 
strongest nebular contamination. Therefore, a larger region of 
5 A was removed from its core. 

The lines from the FLAMES data that will be fitted are the 
hydrogen Balmer lines Ha Hy and H<5; the He i singlet line at 
4387 A; the He i triplet lines at 4026, 4471 and 4713 A, where 
the first line is actually a blend with He n; and the He n lines at 
4200, 4541 and 4686 A. In the formal test we also fit this set 
of lines. Table|2]lists, for the above set of lines for the different 
data sets the final fit parameters obtained by evolving a popu- 
lation of 72 fastwind models over a course of 200 generations. 
Also listed in Tab. |2] are the parameter ranges for each param- 
eter in which the fitting method was allowed to search for the 
best fit. 

For all three tests the automated method was able to find 
the global optimum. There are some differences between the 
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Table 2. Results of the formal tests. Input parameters of the formal test models are given in the "In" column and parameters 
obtained with the automated fitting method by fitting synthetic data created from these models are listed in the "Out" column. 
The results were obtained by evolving a population of 72 fastwind models over a course of 200 generations. 
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0.98 


Voo [kms -1 ] 


3000 






2200 






2000 







final fit parameters and the input parameters. However, these 
can be explained by the low quality of the data simulated and 
the fact that the sensitivity of some parameters becomes less in 
certain parts of parameter space. The latter explains the differ- 
ence between the wind parameters found in case of data set C. 
As the wind of this object is very weak (order 10~ 8 Moyr -1 ) 
very little information about it is available in the spectrum. 
Consequently, the fit parameters describing the wind (J3, M) are 
relatively poorly constrained, with the errors in M becoming as 
large as the actual M value. 

The low S/N also decreases the sensitivity of the fit parame- 
ters. For decreasing S/N the global optimum in parameter space 
becomes shallower. Consequently, the error on the fit parame- 
ters incre ases, as t hese are a measure of the width of this opti- 
mum (see Paper I). This explains the differences in the micro- 
turbulent velocities (~45%), as well as the mass loss rate found 
for data set B. 

To locate the best fit, respectively, approximately 100, 50 
and 60 generations were needed in case of data set A, B and 
C. To obtain robust results we adopt 150 generations as the 
minimum value in fitting our programme stars. 

3.3. Nebular emission 

The subtraction of nebular emission features in multi-fibre data 
can be problematic. The use of a combined sky-spectrum, for 
instance, does not always result in a complete removal of the 
nebular component and in other cases may result in an over- 
subtraction of nebular features. To best cope with these poten- 
tial problems in the automated fitting method we only consider 
the wings of the line profiles, ignoring any core nebular con- 
tamination. As in most cases the co re nebular feature is well 
resolved (also see Evans et al. 120051) and as the removal of the 
line cores does not hamper an accurate determination of the fit 



parameters (see Sect. 13. 2\ . this approach seems justified. It is 
much more difficult to assess the extent to which the determi- 
nation of the fundamental parameters might be influenced by 
residual nebular contamination or subtraction effects in the line 
wings. Here we can only perform limited tests. To assess the 
effect of too much sky-subtraction we looked in more detail to 
the non-sky-subtracted data for NGC 346-010 and NGC 346- 
077. 

To determine if the derived mass loss rates are affected by 
over-subtraction of nebular features in the Ha Balmer line, 
we refitted the non-sky-subtracted spectrum of the 07 giant 
NGC 346-010. From the targets observed with FLAMES this 
object has the smallest mass loss rate that we could deter- 
mine accurately, i.e. with error estimates smaller than 0.2 dex. 
Therefore, if incorrect sky-subtractions would be an issue, 
the mass loss rate estimate of this object would be affected 
the most. In comparison with the fit parameters determined 
from the sky-subtracted spectrum, no significant differences 
are found for any of the parameters obtained from the non- 
sky-subtracted spectrum. In more detail, the mass loss derived 
from the non-sky-subtracted spectrum is found to be larger by 
the small amount of ~0.07 dex, and can be attributed to the 
slightly larger log g (-0.08 dex). As no UV spectrum is avail- 
able for this object, was scaled with u esc , resulting in an ap- 
proximately ten percent higher terminal flow velocity. As M is 
connected to Voo through the continuity equation, this explains 
the 0.07 dex increase in the mass loss rate. 

The second test we performed was on the spectrum of the 
09 dwarf NGC 346-077. This relatively faint object suffers 
from quite severe nebular contamination in its line profiles. 
We find that the fit parameters obtained from the non-sky- 
subtracted spectrum compare well with the values determined 
from the sky-subtracted spectrum. Within the error bars the two 
parameter sets, again, are in agreement. Only a relatively large 
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difference is found for the microturbulent velocity. This param- 
eter was found to be reduced by ~14 kms -1 in the fit of the 
non-sky-subtracted spectrum. However, we do not attribute this 
relatively large change to issues due to nebular contamination. 
Instead, as the formal tests have shown, an accurate determina- 
tion of this parameter is notoriously difficult for low signal to 
noise spectra (see also Sect. 14. 5> . 

Now, what if our combined sky-spectrum underestimates 
the real sky background in the line wings? To assess this, ide- 
ally one would like to compare the current fitting results to 
results obtained from data with a local sky-subtraction. For 
NGC 346-001 we had the opport unity to perform such a test, 
as also the spectrum analysed by ICrowther et alJ (|2002) was 
available. A fit of this spectrum resulted in parameters nearly 
identical to those obtained from the VLT-FLAMES spectrum. 
In particular the wind parameters /? and M, which are expected 
to be the most sensitive to nebular contamination, were found 
to agree within, respectively, 10 percent and 0.03 dex. In itself, 
this agreement is reassuring, though, we note tha t NGC 346- 
001 i s located away from the cluster centre (see lEvans etal] 
2006). Consequently, other objects might still suffer more from 
nebular contamination. With respect to the mass loss rate de- 
terminations, we note that from the objects in NGC 346 only 
four have a reliable mass loss rate. Out of these only one object 
(NGC 346-033) lies close to the core of the cluster. Therefore, 
its mass loss may be more uncertain than is suggested by 
the formal errors (also see Sect. 14. 6> . The other three objects 
(NGC 346-001, -010 and -012) lie at relatively large distances 
from the core, where we anticipate the background contribution 
to be relatively small. 

3.4. Error estimates 

We define error estimates for the fit parameters by estimat- 
ing the width of the optimum in parameter space associated 
with the global optimum. This width defines, as was argued in 
Paper I, the region in parameters space which contains mod- 
els with comparable fit quality. Consequently, by determining 
the maximum variation of the par ameters within this group of 
models, the error is estimated (see Paper!!). 

Table |4] contains the optimum width based error estimates 
of the fit parameters for each analysed object. Based on these 
estimates the errors on the derived parameters (/?+, L*, M s and 
M ev ) in this table were calculated using the same approach as 
was used in IPaoer 1 The single difference is the adopted un- 
certainty in the absolute visual magnitude. Here we adopt an 
uncertainty of 0.14™. This value is equal to the sum of the sta- 
tistical and systemat ic error in the deter mination of the SMC 
distance modules bv lHarries et al.l J2003h . The method used to 
determine the uncertainty in go is explained in Sect. 14.21 

4. Fundamental parameters 

In this section we will discuss the stellar properties of the in- 
vestigated sample. Table [3] lists the values determined for the 
seven free parameters as well as quantities derived from these. 
Error estimates on the parameters are given in Tab. |4] The fits 
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spectral type 

Fig. 1. Effective temperatures as a function of spectral type for 
SMC objects studied in this paper. Dwarfs, giants and super- 
giants are denoted by, respectively, solid circles, triangles and 
squares. Shown with a solid line is the average T e ff relation for 
the dwarfs with spectral type later than 04. SMC dwarfs are 
found to be typically 3 to 4 kK hotter t han the 7V.ff cal i bration 
for Galactic O-type dwarfs according to Martin s et alJ (12005a. 
dashed lin e). The dotted line co rresponds to the SMC T e ff scale 
derived by M assev et alJ d2005l) for luminosity class V and III 
objects. 

of the spectra together with comments on the individual objects 
are presented in the appendix. 

4.1. Effective temperatures 

The cumulative opacity of all spectral lines, referred to as line 
blanketing, has a strong effect on both the structure and the 
emergent spectrum of hot star atmospheres. It has been shown 
by several authors that line blanketing changes the relation be- 
tween spectral type (set by the ionisation balance of mainly 
helium and/or silicon) and effective temperature (related to the 
gas temperature of the line forming layers). As line blanket- 
ing enhances the diffuse radiation field, because lines "trap" 
the photons and introduce additional back scatte ring (see e.g. 
ISchaerer & de Koterlll997t iRepolust et aDl2004 . models that 
account for it can su ffice with a lower temperature to match a 
given spectral type jde Koter et al" Il99gt iMartins etal] [2002: 
ICrowther et"aDl2002HHerrero et al.ll2002h . As iron group lines 
dominate the line blanketing, the spectral type vs. T e ff relation 
is expected to depend on metallicity, i.e. SMC stars of given 
spectral type will have a higher effective temperature compared 
to their galactic counterparts. 

Figure^shows the distribution of the effective temperature 
as a function of O and early-B spectral sub-type for the in- 
vestigated sample. The different luminosity classes are denoted 
using circles for the dwarfs and subgiants, triangles for the gi- 
ants, and squares for the bright giants and supergiants. We now 
concentrate on dwarf stars only, which account for the vast ma- 
jority in our sample. For these objects the figure shows a well 
defined relation, the mean of which is represented by the solid 
line. At the earliest spectral types (04-06) the relation seems 
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Table 3. Fundamental parameters determined using GA optimised spectral fits, with r e ff in kK, log g and log g c in cm s~ 2 , R+ in 
Rq, in Lq, w tur b and u r sin/ in kms -1 , M in M0yr _1 , M s and M ev in Mq and Qo in number of photons per second. Results 
were obtained using a population of 72 fastwind models evolved over a minimum of 150 generations. Gravities corrected for 
centrifugal accele ration (logg c ) were used to calculate the spectroscopic masses (M s ). Evolutionary masses (M ev ) were derived 
from the tracks of lCharbonnel et alJ Jl993h . Note that "NGC" is omitted from the identifications of the FLAMES targets. 
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A value of ji = 0.80 corresponds to an assumed fixed value 



to flatten out. However, the small number of objects analysed in 
this range make this part of the diagram uncertain. Moreover, 
two of the objects at spectral type 04 and 05 suffer from strong 
nebular contamination, complicating the determination of their 
effective temperature (see the appendix). 

Shown with a dashed line in Fig. ^ is the T^g "observa- 
tional" calibratio n for G alactic O-type dwarf stars as derived 
by iMartins et alJ ( 12005 al) . A clear offset is apparent between 



this relation and the average calibration for the SMC O dwarfs. 
This offset is approximately 3.3 kK for the latest types up to 
approximately 4.4 kK for spectral type 05. These differences 
in T e ff correspond to a shift in spectral sub-type of about 2 for 
the late type objects a nd ~1.5 sub-typ es for the earliest types . 
Moki em et all |2004), Using CMFGEN faillier&Milleri ri998) 
and covering the metallicity range from 2 times Solar to 1/10 
Solar, find typical shifts of one sub-type, though we should add 
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Table 4. Optimum width based error estimates for the seven fit parameters. The ND entries correspond to error in i> tur b that reach 
up to the maximum allowed value of u tul b an d, therefore, are formally not defined. See text for details on the calculation of the 
uncertainties in the derived parameters. Units: r e tj in kK, log g and log g c in cm s~ 2 , R+ in Rq, L+ in Lq, v tm b and v r sin i in km s _I , 
M in Mq yr _1 , M s and M ev in Mq and Qo in number of photons per second. Note that "NGC" is omitted from the identifications 
of the FLAMES targets. 
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that this comparison is not as extensive as the one presented 
here. 

Figure [2 also shows a clear separation between objects 
of different luminosity class. Compared to dwarfs, the giants, 
bright giants, and supergiants systematically have lower effec- 
tive temperatures. The reason for this separation is twofold. 
First, the latter group of objects represent more evolved evolu- 
tionary phases. Thei r lower gravities resul t in an increased he- 
lium ionisation (e.g. Kudritz ki et alJll9 83). reducing the effec- 



tive temperature associated with a given spectral type. Second, 
these objects are expected to have stronger stellar winds. This 
induces an increased line blanketing effect, further reducing the 
r e ft for a given spectral type. 

Mass ev et alJ 12005) also report a Sp.Type(r e +r) calibration 
for SMC stars. As can be seen in Fig.^their relation essentially 
agrees with ours at spectral types earlier than 08. However, at 
later types the ir results suggest a rather sharp turn towards the 
Mart ins et alJ relation for Galactic stars which is not observed 
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Fig. 2. Number of Lyman continuum ionising photons Qo as 
a function of spectral type. Different luminosity classes are 
shown using circles, triangles and squares for, respectively, IV- 
V, III and I-II class objects. Grey circles correspond to stars lo- 
cated on or left of the ZAMS. Indicate d using a dashed line i s 
the Galactic calibration for dwarfs from M artins et alJ J20 05a). 

in our sample. The reason for this apparent discrepa ncy is that 
our calibration employs dwarfs, whilst Mas sev et alJ had to rely 
upon giant stars at the latest O subtypes. 

4.2. Ionising fluxes 

The ionising fluxes of massive stars are important quantities 
that are used in the st udy of, for instance, H n regions and star- 
burst galaxies (e.g. IVaccalll994l) . The parameter that is used 
to characterise the ionising output is the number of photons 
present in the Lyman continuum, Qo- It is defined as 

Qo = 4^ C^dA, (1) 
Jo he 

where and Ao are, respectively, the stellar flux in 
erg s _1 crrT 2 A -1 and the limiting wavelength for photons able 
to ionise hydrogen, i.e. 912 A. 

Tables|3]and|4]list the number of ionising photons for the in- 
dividual objects and the associated error estimates. These errors 
are dominated by the uncertainty in stellar radius and effective 
temperature. To calculate the error estimates we adopted an un- 
certainty of 0.03 dex in which is dominated by the estimate 
for AMy. Th e error introduced by AT e ff was estimated from 
Fig. 16 from lMartins et all J2005al) . which shows the number 
of ionising Lyman continuum photons as a function of effec- 
tive temperature for the line blanketed stellar atmosphere codes 
cmfgen, wm-basic and tlusty. The difference between the pre- 
dictions of these three codes is relatively small and, conse- 
quently, we estimated the uncertainty introduced in Qo to be 
0.09 dex/kK for r eff > 36 kK and 0. 1 8 dex/kK for r eff < 36 kK. 

In Fig. 13 the distribution of Qo as a function of spectral 
type for our programme stars is presented. Different luminos- 
ity classes are indicated using circles, triangles and squares for 
respectively, IV- V, III and I-II class objects. On average the 
more evolved stars are found to produce more ionising pho- 



tons for a given spectral type because of their larger radii, 
hence higher luminosities. Also shown in this figure as a dashed 
line is the "observational" Qo calibration for Galactic dwarfs 
from lMartins et al.1j2005al) . In general we find that the ionising 
fluxes of the SMC dwarfs are in good agreement with this cal- 
ibration (also see iMokiem et alJl2004t) . The most pronounced 
differences are for the 04 and 07 stars. With respect to the 
earliest spectral type this can be explained by nebular contami- 
nation hampering the T e ff determination (see previous section). 
The low ionising fluxes of the two 07 dwarfs, however, can- 
not be explained in this way, as they are found to follow the 
r e ff trend in Fig.^ Instead, we believe that the discrepancy is 
related to their location in the HR-diagram, i.e. evolutionary 
phase. The two dwarfs are found to locate a position on or to 
the left of the ZAMS (see Sect. 17121. As a result of this they are 
less luminous compared to objects of the same spectral type 
that are located to the right of the ZAMS. Consequently, the 
total number of ionising photons produced by these objects is 
smaller than the average associated with their spectral type. 

Apart from the two dwarfs at spectral type 07 three addi- 
tional stars were found to lie on or to the left of the ZAMS (see 
Sect. 17. 2\ . As can be seen in Fig. 13 where we have highlighted 
all ZAMS stars using grey circles, these stars produce on aver- 
age less ionising photons. Note that two dwarfs at 08 and 08.5 
also seem to lie below the average. The 08.5 star has a below 
average T e ff for its spectral type (see Fig.0. This explains its 
somewhat peculiar Qo behaviour. The 08 star below the av- 
erage is NGC 346-050. It also has a temperature that is lower 
than the average for its spectral type, though not to the same 
extent as NGC 330-052. Interestingly, this object lies closest to 
the ZAMS (see Fig. [O} of all non-ZAMS stars, and seems to 
behave in terms of Qo in a similar manner as the ZAMS ob- 
jects. We conclude that the ZAMS stars for given spectral type 
have ~0.4 dex lower Lyman continuum photons. 

4.3. Gravities 

In Fig. [3] we present the distribution of SMC objects in the log 
T e ff - log g c plane. The surface gravity corrected for centrifugal 
acceleration ( logq r ) was c a lculat ed ac cording to the metho d 
discussed by Herre ro et a ill 9921) and lRenolustetal] {2004). 
Different luminosity classes are indicated using circles, trian- 
gles, and squares for type IV- V, III, and I-II objects, respec- 
tively. With the exception of one object, the dwarf type objects 
are all above \ogg c — 3.9. The object with a lower gravity, 
NGC 346-026 is the only subgiant. The luminosity class I-III 
objects occupy a strip below the dwarfs, reflecting the evolu- 
tionary path of hot massive stars in this diagram. Note that no 
separation between luminosity class III and I-II objects is visi- 
ble in this secular behaviour. 

A comparison of the masses based on the spectroscopi- 
cally determined surface gravities and those derived from pre- 
dictions of massive star evolution is presented in Fig. |4] The 
different luminosity classes are distinguished using the same 
symbols as in Fig. [5] To deter mine the ev o lution ary masses 
evolutionary tracks from Char bonnel et alJ dl993h . for Z = 
0.2 Zq were used. The errors in evolutionary mass reflect the 
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Fig. 3. Distribution of the analysed SMC objects in the log T e s 
- log g c plane. Indicated with different symbols are luminosity 
class IV-V (circles), III (triangles) and I-II (squares). With ex- 
ception of the subgiant NGC 346-026, the first group of objects 
are clearly separated from the more evolved objects. These 
more evolved objects seem to follow a trend of increasing sur- 
face gravity with increasing effective temperature. 



mass interval allowed within the error box spanned by the 
stellar luminosity and effective temperature. As the tracks of 
ICharbonnel et al.l do not account for the effects of rotation, 
this source of error is not included. Predictions accounting for 
v Iot show complicated tracks including loops during the sec- 
ular redward evolution. Therefore, one can no longer assign 
an unambiguo us M(L, T^). Still , asses sing the impact of rota- 
tion us ing the lMaeder & Mevnetl JboOl ) and lMevnet & Maedell 
(2005) computations that adopt an initial rotational velocity 
Vj = 300kms _1 shows that the error in the evolutionary mass 
will not increase by more than ~13 percent. The errors in the 
spectroscopic mass are much larger than those in M ev , and pri- 
marily reflect the error in gravity. 

Inspection of Fig. |4]reveals no convincing systematic dis- 
crepancy between the spectroscopic and evolutionary mass, 
even though some objects do not agree within their standard de- 
viation with the one-to-one relation. At the low mass end there 
appears to be some tendency for the spectroscopic masses to be 
less than those from evolutionary tracks. This behaviour is sim- 
ilar to t he "mass discrepancy pro blem" reported and dis cussed 
by e.g. iHerrero et all d2002l) and iRepolust etaD d2004 . Note 
that most of this classical problem has been resolved, i.e. it has 
been attri buted to limitations of the stellar atmospher e model s 
jHerrerolll993l) and biases in the fitting process (see lPaperll) . 
With respect to the stars at the low mass end, we note that star 
NGC 346-107 occupies a location in the HR-diagram left of the 
ZAMS (see Fig.|6}. As the evolutionary status of this object is 
formally not defined, the mass that is given is based on an ex- 
trapolation of the tracks. This could lead to an erroneous value 
of M ev . Interestingly, for the other stars at the low mass end the 
helium abundances listed in Tab. [3] seem to correlate with the 
mass discrepancy. We will investigate this in detail in the next 
section. Also note that Massev et al. I ll2005l) identify a mass dis- 
crepancy in a sample of Magellanic Cloud stars for objects with 
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Fig. 4. Comparison of spectroscopically determined masses 
with masses derived from the evolutionary tracks of 
CharbonniriLal] l<1993l) for Z = 0.2 Zq. Denoted by circles, 
triangles and squares are, respectively, luminosity classes IV- 
V, III and TIL The dashed line corresponds to the one-to-one 
relation between the two mass scales. With exception of the 
two most massive objects, no systematic discrepancy between 
the spectroscopic and evolutionary masses is observed. 



T e g > 45 kK, which they attribute to a possible underestimation 
of log g by the model atmospheres. Unfortunately, our sample 
does not contain enough early-type O stars to corroborate the 
findings of these authors. 

For the two brightest objects, AzV 14 and AzV 26, the 
spectroscopic masses are much larger than the implied evo- 
lutionary masses. The profile fits, in particular those of the 
gravity sensitive hydrogen Balmer lines, (see the appendix) are 
good. So, it is not likely that the reason for the discrepancy 
is an overestimation of the spectroscopically determined grav- 
ity. One could speculate about a possible binary nature of both 
stars, as the spectroscopic mass is more sensitive to changes in 
luminosity. An indication of binarity may be the conspicuously 
high luminosities. Howev er, using spectroscop ic and spectral 
morphological arguments, Massev et alJ J20 04) tentatively rule 
out a composite explanation for both objects. In contrast to this 
we note that the luminosity and derived mass-loss rate for these 
two stars (see Sect. I4.6I > imply a position in the modified wind 
momentum vs. luminosity diagram which is well below what 
is expected from theory, i.e. the mass loss for these two stars 
would be in better agreement with predictions if their luminos- 
ity would be lower. 

4.4. Helium abundances 

The automated method also treats the helium abundance as a 
continuous free parameter. Usually, spectroscopic analyses as- 
sume an initial Solar value for Yn e , whi ch is only mod ified 
when no satisfying fit can be obtained (e. g. lHerrero et alJ2 002: 
IRepolust et alJ2004l:lMassev et al.l2004 . Because of the auto- 
mated treatment and the extent of our sample we can, for the 
first time rigorously, investigate possible correlations between 
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Fig. 5. Helium abundances as a function of surface gravity. 
Luminosity classes IV- V, III and I-II are denoted using cir- 
cles, triangles and squares, respectively. Open symbols corre- 
spond to fast rotators (u r sin/ > 250 km s -1 ). The average he- 
lium abundance of the investigated sample is 0.13. Shown with 
a dashed line is a measure for the initial Fee, which is taken to 
be the average of the helium abundances of the dwarfs with a 
Ine smaller than the sample average. Compared to this "initial" 
abundance of 0.09 an increase in Yn e is visible for decreasing 
surface gravity. 




Fig. 6. HR-diagram of the analysed SMC sample. Different 
luminosity classes are distinguished using identical symbols 
as used in Fig. |5] Open symbols correspond to objects 
with a helium abundance of at l east 0.12. Over-plo t ted in 
grey are the evolutiona r y track s of iMaeder & Mevnell (12001 ) 
and iMevnet & Maederl (120051) with Z = 0.2 Z and v r = 
300 km s -1 . The ZAMS corresponding to these tracks is shown 
as a black solid line. The grey area in the top part of the diagram 
corresponds to the region in which the evolutionary models ex- 
hibit an enhancement of the surface helium abundance by ten 
percent or more. 



the surface helium abundance and other fundamental parame- 
ters. 

In Fig. |3 we show Yu e as a function of surface gravity 
log g c . The horizontal dashed line represents the initial helium 
abundance of the SMC stars investigated. It corresponds to the 
average of the helium abundances of the dwarf objects with 
Fne smaller than the total sample average. Using this "initial" 
abundance of 0.09^'^ as a reference, the overall trend is that 
the average Yu e increases for decreasing surface gravity. This 
is consistent with the standard picture that more evolved ob- 
jects may have their atmospheres enriched with primary he- 
lium. Interestingly, however, one may immediately spot two 
deviating objects from the overall trend. First, the supergiant 
NGC 346-012 has a helium abundance lower than the "initial" 
value. This object has \ogg c = 3.35 and Yu e = 0.07. The rea- 
son for the low surface helium abundance is unclear and we 
will exclude this star from the remainder of this discussion. 

Second, some of the unevolved objects have enhanced 
helium abundances. It can therefore be suspected that more 
parameters are involved in controlling the enrichment dis- 
play ed in Fig. El One such parameter could be stellar rotation. 
Meynet & Maeder (2000), for instance, predicted that exten- 
sive mixing in fast rotators could result in significant surface 
helium enhancement relatively early in the evolution. To probe 
this possibility we have highlighted the fast rotators, defined as 
having u r sin/ > 250kms~', in Fig. fusing open symbols. In 
case of all four fast rotators we see that their helium abundances 
are enhanced with respect to the "initial" value, suggesting that 
the helium enhancement may be (partly) related to fast rota- 
tion. This may imply that two other dwarf objects with a clear 



helium enhancement (Yn e » 0.16) but with low projected rota- 
tional velocity are in fact fast rotators seen pole on. 

How well do our spectroscopically derived helium abun- 
dances compare to those predicted by evolutionary models? 
This is a fundamental but complicated question. Important to 
realise is that one of the effects of rotation is that it introduces a 
wide bifurcation in the evolutionary tracks. Stars rotating above 
a threshold v t o f about 30-50 percent of break-up (depending 
on initial mass; lYoon & LangeJ2005l) follow tracks which are 
essentially those of chemically homogeneous evolution. Such 
tracks do not evolve towards the red in the HRD, but evolve 
bluewards (from the zero age main seque nce) and upwa rds un- 
til the star enters the Wolf-Rayet phase (Maede3 ll987l) . Stars 
rotating below this critical value evolve along tracks which are 
about similar to non-rotating ones, though at an earlier age the 
surface abundances of helium (and carbon, nitrogen and oxy- 
gen) will be affected by rotation induced mixing. 

Let us first compare our results with tracks for rotational 
velocities below this critical value. In the HR-diagram shown 
in Fig. the g rey lines corres pond to evolutionary prediction s 
by |Maeder & MevneJ J200lll and IMevnet & Maeder ( 2005), 
which were calculated for Z = 0.2 Zq and v rot = 300kms _1 . 
Objects for which a spectroscopic helium abundance of at least 
Yhss = 0.12 was found are denoted using open symbols. Note 
the location of four helium rich dwarf stars close to or even on 
the ZAMS. These will be discussed in more detail in Sect. 
The grey area in the figure corresponds to the region in which 
the evolutionary models predict a surface helium enhancement 
of at least ten percent. 
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Fig. 7. The mass discrepancy (see also Fig.|4} as a function of 
surface helium abunda nce. To compute th e discrepancy non- 
rotating tracks of Charbonnel et al. (1993) are used. Circles, 
triangles, and squares denote dwarfs, giants, and supergiants, 
respectively. The open circles indicate stars on or left of the 
ZAMS. The bulk of the stars has the initial helium abundance 
^He = 0.09 where the discrepancy seems random and of the 
order of the error. For helium enriched stars the discrepancy al- 
ways favours a higher M ev , qualitatively consistent with models 
accounting for mixing increasing the M/L-ratio and bringing 
helium to the surface. A quantitative prediction (short dashed 
lines) based on ZAMS stars of 20 (upper curve) and 30 Mq 
(lower curve) varying helium abundance shows, however, a 
more modest effect than implied by our findings. 

The regime in which the evolved stars showing significant 
helium surface enrichment reside roughly coincides with the 
location of the grey area. Not all evolved objects show evidence 
of Yn e enrichment. Many exhibit an abundance about equal to 
the initial value. None of these stars are fast rotators (they all 
have u r smi < 140kms _1 ). This can easily be explaine d us- 
ing tracks for non-rotating stars (Mevnet & Maeder 2000). The 
fair consistency between observed and predicted helium abun- 
dance in evolved objects is reassuring, however, a more detailed 
comparison requires the availability of tracks for several more 
values of u lot . 

We mentioned above that for modest rotation the tracks 
do not differ greatly from non-rotating ones. This is not com- 
pletely correct as rotatio n tend s to make the star somewhat 
more luminous. Langerj i 19921) proposed that this mig ht ex- 
plain t he mass discrepancy problem identified bv lHerrero et alJ 
Although recent analyses no longer suggest a convinc- 
ing systematic discrepancy (see e.g. Sect. l4.3t . it is still interest- 
ing to look at this idea in some more detail. Langer connects the 
apparent mass problem to the helium abundance by showing 
that the M/L-ratio is a monotonically decreasing function for 
increasing helium enrichment. Consequently, if mixing brings 
primary helium to the surface one expects the star to be overlu- 
minous, leading to an overestimate of the evolutionary mass if 
non-rotating tracks are adopted. If indeed this is the case, then 
the scatter around the one-to-one relation in Fig.0]might reveal 
a trend when plotted against Yft e . 



The result of this exercise is shown in Fig.0 On the verti- 
cal axis a measure of the mass discrepancy is given. Note that 
the mass difference is plotted relative to the mean of the evolu- 
tionary and spectroscopic mass to ensure that positive and neg- 
ative discrepancies are shown on the same linear scale. For the 
evolutionary masses, tracks that do not account for rotation are 
used. The circles, triangles, and squares denote dwarfs, giants, 
and supergiants, respectively. The open circles indicate stars 
on or left of the ZAMS. At the "initial" helium abundance of 
iHe = 0.09 + 0.01 - where most of the stars reside and dwarfs 
dominate - the scatter around the M ev = M s relation appears 
random. This essentially reflects that there is no systematic 
mass discrepancy. At 0.13 < Tne $0.19 the scatter is not ran- 
dom, as all objects show a positive mass discrepancy. In prin- 
ciple this is qualitatively consistent with the above described 
idea. However, is it also quantitatively consistent? To assess 
this we have computed the mass discrepancy for ZAMS stars 
with a variable helium abundance. This should reflect the max- 
imum effect of rotation, i.e. such effective mixing that it leads 
to chemically homogeneous evolution. The results for stellar 
masses of 20 and 30 Mq - typical for the bulk of our sample - 
are shown (short dashed lines). These predictions clearly show 
a more modest mass discrepancy, though the error bars on the 
mass discrepancy for the programme stars do reach these pre- 
dictions. We conclude that stars with an enriched helium sur- 
face abundance tend to show a systematic mass discrepancy 
that is qualitatively consistent with predictions of chemically 
homogeneous evolution. We finally note that the supergiants in 
our sample, which can be explained using evolutionary models 
including rotation, show the best agreement between M ev and 
M s (see also Fig.0}. 



4.5. Microturbulence 

Similar to the helium abundance the fact that the microturbu- 
lent velocity is treated as a free parameter allows us for the 
first time to investigate possible correlations for this parameter. 
However, in contrast to the former parameter, this investigation 
did not yield any clear relation between i> tm b and any other pa- 
rameter. For each parameter given in Tab.|3]a comparison with 
the microturbulent velocity basically results in a scatter dia- 
gram. This null result is similar to the findings in IPaoer 3 and 
reflects the uncertainty with which u tul -b can be determined from 
the hydrogen and helium spectrum. Apparently the line profiles 
are not very sensitive to this parameter. 

We should also consider the fact that the error estimates 
determined for i> tur b, as given in Tab. |4] are on average consid- 
erably large. As a result of this, one could argue that for many 
objects it simply was not possible to accurately determine i> tur b. 
Consequently, it is not possible to find any correlation when 
the total sample is considered. To avoid this potential problem 
we also investigated possible correlations using a subset of 14 
objects for which the microturbulent velocity was determined 
relatively well. The selection criterion for this subset was that 
the Uturb error bars should be well confined within the search 
domain, spanning the range of up to 20 km s _1 . The compar- 
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Fig. 8. Modified wind momentum (D mom ) in units of gems' 2 Rq as a function of luminosity for the analysed SMC sample. 
Luminosity classes are distinguished using identical symbols as in Fig. [5] Indicated using grey inverted triangles are upper limits. 
At log L+/Lq ~ 5.0 the open symbol indicates the location NGC 346-033 would have if its true v m would be lower by a factor 
two compared to the adopted value. Using open squares modified wind momenta corrected for clumping are shown for the 
two supergiants with Ho- emission profil es. The upper and lowe r dashed lines correspond to the theoretical wind-momentum 
luminosity relations (WLR) as derived bv lVink et allfeOOOl koOl ) for, respectively, Galactic and SMC metallicity. Shown with a 
dotted line is the empirical WLR constructed for objects with log L+/Lq > 5.4. 



ison within the subset again did not result in a correlation for 
^turb with any of the other parameters. 

4.6. Wind parameters 

Our sample is dominated by late O-type dwarf stars, which are 
expected to have relatively weak winds. These are so weak that 
they challenge the sensitivity of Ho- as a mass-loss diagnostic. 
With nebular contamination as an additional complicating fac- 
tor, we could not derive reliable M for all objects. For nineteen 
stars (see Tab.|5J we can only derive upper limits, i.e. the down- 
ward error bars on our fit extend to the lower limit of the M 
regime in which the automated method was allowed to search 
for a solution. It happens to be so that these can be identified 
by an error bar - log AM > 1.0 dex (see Tab. [3}- For many 
stars we could also not determine the acceleration behaviour of 
the wind, expressed by the exponent £1 of the velocity law. For 
these objects we adopted = 0.8 consistent with theoretical 
expectations (Paul drach et alJl986l) . 

Despite the large number of upper limits, it is still pos- 
sible to quantitatively investigate the SMC stellar winds. We 
do this by studying the distribution of the analysed objects in 
the so-called modified stellar wind momentum vs. luminosity 
diagram. The modified wind-momentum, which is defined as 
D mom = Mv m R]J 2 , is predicted to behave as a power-law as a 



functi on of stellar luminosity JKudritzki et al]ll995l IPuls et all 
Il996h 

log D mom = x log (L+ /Lq) + log D . (2) 

In the above equation x corresponds to the inverse of the slope 
of the line-strength distribution function corrected for ionisa- 
tion effects (IPuls et alJ l2000) and D is a measure for the ef- 
fective number of lines contributing to the acceleration of the 
wind. 

In Fig.|S]we present the distribution of the modified wind 
momenta of our sample. The upper dashed curve is the theoret- 
ical prediction for a Galactic metal content; the lower dashed 
curve is that for an SMC metallicity (iVink e t al. 2001; see also 
below), which is predicted to be shifted downward by 0.57 dex 
with respect to the Galactic relation. Before confronting the- 
ory with observations, we first discuss a few individual objects. 
The D mom upper limit at logL+/LQ » 5.9 corresponds to the 
object AzV 14. As was argued in Sect. 14. 31 this object might be 
over-luminous because of a binary nature. 

The 08 V star NGC 346-033 (at log L*/L * 5.0) is 
positioned far above the SMC prediction. We suspect this 
is connected to the anomalously high terminal velocity of 
4100 kms _I , which results from a scaling with u esc as no di- 
rect UV measurement is available. If this lie would be overes- 
timated by a factor of two (for an 08 V star one would expect 
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Uoo ~ 1900 km s"\ cf. lKudrit/ki &Piildl200fih the D mom would 
be reduced by approximately a factor of six. This is so because 
D mom scales directly with v m and indirectly with v^ 12 through 
the invariant wi nd-strength parameter (Q oc M/v^^ 2 , e.g. see 
IPuls et alJ2005l) . The open circle shows the effect of such a de- 
crease in terminal velocity. We decided to exclude this object 
from the remainder of this discussion. 

The two O-type supergiants in our sample, AzV 372 and 
NGC 346-00 1, are the only two exhibiting an Ha em ission 
profile. Mark ova et alJ J2004 and lReoolust et alJ J2004 have 
argued that for such stars the mass loss may be overesti- 
mated relative to dwarf stars due to wind clumping effects. 
For dwarf stars the Ho- absorption line is formed relatively 
close to the stellar surface, where clumping may be negli- 
gible. The Ha emission line is typical for supergiants and 
in contrast reflects that the line is formed over a larger vol- 
ume, where - the y propose - clumpin g has set in. Indeed, 
for NGC 346-001 ICrowther et all J2002I) present evidence for 
wind clumpin g based on th e analysis of the UV phosphorus 
lines. iRepolust et all J20 04) derive a correction for this clump- 
ing by multiplying the mass loss by a factor 0.44, which we 
have applied for the two objects. The corrected modified wind- 
momenta are shown in Fig. fusing open squares. 

For SMC objects, Fig. [8] represents the best populated 
modified wind-momentum diagram presented so far (12 mass 
loss determinations and 19 upper limits). In particular at 
log L+ /Lp, > 5.4 we find excellent agreement with the 
IVink et alJ predictions, establishing that at these luminosities 
the winds of SMC stars are weaker, in accord to theoretical ex- 
pectations. At log L+/Lq < 5.2 the situation is less clear as 
for many stars we could only set upper limits and for those for 
which M could be derived the error bars are large. For the weak 
wind regime we therefore cannot draw firm conclusions. 

We have constructed an empirical WLR by fitting a lin- 
ear function to the objects in Fig. [8] with log L+/Lq > 5.4, 
while accounting for the symmetric errors in luminosity and 
the asymmetric errors in D mom . Using the clumping corrected 
£>mom values for AzV 372 and NGC 346-001 this results in the 
following relation 

logD mom = (1.69 ± 0.34) log (LJLq) + (18.87 ± 2.00) , (3) 

which is shown in Fig. [8] as a dotted line. Within the error 
bars the fit pa rameters agree with the theoretical calculations 
of lVink et alJ feOOlb . who predict x = 1.83 and log D a = 18.11 
for SMC metallicity. To illustrate the agreement we compare 
the predicted and fitted wind momenta at the boundaries of 
the fitting range. For log L+/Lq = 5.4 and log L+/Lq = 6.2 
the differences are -0.02 dex and 0.08 dex, respectively. Given 
the fact that the typical uncertainty in D mom is of the order of 
0.2 dex, we find the agreement excellent. We note that the weak 
wind regime is discussed in more detail in Sect. [6] 

4.7. Projected rotational velocities 

The projected rotational velocities as determined with the au- 
tomated fitting method are listed in Tab. [3] The associated un- 
certainties range from approximately 10 kms -1 for the slow 
rotators up to about 40 kms -1 for the fast rotators, indicating 



that we were able to accurately determine u T sin i with our self 
consistent method. Note that the low u T sin i of NGC 346-03 1 
and NGC 346-051 are close to the effective resolution of the 
FLAMES observation (~15 kms 1 ), and should be interpreted 
as upper limits. 

If a star rotates at more than approximately 80% of critical 
rotation the assumption of spherical symmetry may break down 
(e.g. Maeder & Mevnet 2000). Checking the ratio v r sin i/v c „ t 
for all our target stars, we find a maximum of about ~0.65. We 
therefore do not expect any significant deviations from spheric- 
ity, nor uncertainties in the derived parameters due to rotational 
effects. 

4.7.1. Rotation vs. macroturbulence 

In our calculations we consider the possibility of small scale 
velocity fields. These thermal and microturbulent motions have 
a coherence length that is small compared to the line form- 
ing region. We do not include the possibility of motions that 
have a coherence length that is comparable to or is larger than 
the line forming zone. Such motions are denoted with the term 
macroturbulence (v mac ). If a large macroturbulent velocity field 
is present, it may be expected that the projected rotational ve- 
locity that we determine from the hydrogen and helium profiles 
is overestimated. 

To assess the importance of macroturbulence we have 
also determined v r sini by means of a Fourier technique 
dSimon-Dfaz et alJl2006h using weak metal lines as an inde- 
pendent diagnostics. The Fourier technique allows to discrim- 
inate between macroturbulence and rotation, as it is sensi- 
tive to differen ces in the rotational and macroturbulent profiles 
(see e.g. lGravHl978l) . For details on this method we refer to 
Sim on-Diaz et alJ In all but one case we get u r sin i values that 
are consistent to within the error estimates. As in first order 
the Vj sin i and u mac broadening should be added quadratic ally, 
this does not necessarily imply that macroturbulent motions are 
absent. For u r sin i of order 100 to 150 kms~', macroturbulent 
fields with characteristic velocities of up to several tens of kilo- 
metres per second, i.e. comparable to the small scale velocity 
component, remain a possibility. Obviously, rapid rotators may 
show larger u mac components. 

For the 08.5 bright giant AzV 469 the Fourier method re- 
covers a v r sin i that is about 30 percent lower than the 8 1 km s _1 
derived from the hydrogen and helium profiles. This could 
indicate the presence of a macroturbulent velocity field with 
fmac ~ 50kms _I . Significant macroturbulent fields have been 
reported for B supergiants (e.g. lRvans et alJEo02l) . but not for 
early O-type stars. The fact that this star is of relatively late O 
sub-type seems consistent with these findings. Consequently, 
in the u r sin i distribution analysis presented below we adopt the 
projected rotational velocity for AzV 469 as determined using 
the Fourier method. 



4.7.2. Observed i> r sin/ distributions 

For a meaningful comparison of the distribution of the pro- 
jected rotational velocities of the SMC objects with other ob- 
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Fig. 9. Comparison of u r sin ; distributions. Left panel: cumulative distribution functions (cdf) of v T sin i for unevolved objects, i.e. 
luminosity class IV and V, and evolved objects, i.e. luminosity class I, II, III, in the SMC. The group of unevolved SMC objects 
are found to contain relatively both more fast rotators (u r sin i > lOOkms -1 ) as well as more slow rotators (v r sin / < 50km_s__J) 
with r espect to the group of evolved objects. Right panel: cdf of unevolved SMC objects and unevolved Galactic objects JPennvl 
1996). Compared to the Galactic sample the SMC sample contains a larger fraction of intermediately fast rotating stars (120 > 
v T sin i < 190 kms -1 ). 



servations and with theory we use cumulative distribution func- 
tions (cdfs). The cdf describes the distribution of u r sin i by sim- 
ply giving for every observed u r sin i the fraction of objects with 
lower or equal velocities. In Fig.|9]the cdfs of the SMC sample 
are presented and compared to cdfs of Galactic O-type stars. 

The left panel in Fig. |9] compares the cdfs of unevolved 
objects, i.e. luminosity class IV and V, and evolved objects, 
i.e. luminosity class I, II and III, in the SMC. This compari- 
son shows that compared to the evolved objects the group of 
unevolved objects contain relatively more fast rotators. For in- 
stance, approximately ten percent of the evolved objects have a 
n r sin i in excess of 150 km s _1 , whereas approximately 40 per- 
cent of the unevolved objects exhibit velocities larger than this 
y r sin i. Note that about 20 percent of the group of unevolved 
objects is rotating slowly (y r sinj < 50kms _I ), while this is 
only 10 percent for the evolved stars. 

Using the Kolmogorov-Smirnov (K-S) test we have de- 
termined that the probability that the two samples are drawn 
from the same underlying distribution is 23 percent. Therefore, 
the differences between the cdfs of the unevolved and evolved 
SMC stars may be significant and possibly not due to statistical 
fluctuations. The trend that is found here is also seen in Galactic 
O-t ype stars. Us i ng a s imilar approach, iHowarth et alJ dl997l) 
and lPennv et al.l (120 04) also find relatively mor e slow and fast 
rotato r s among the un evolved stars (also see IConti & Eb bets 
Il977t) . IHowarth et al.l ascribe the reduced number of fast ro- 
tators to spin down as a result of an increased radius for the 
evolved stars, as well as to loss of angular momentum through 
the stellar wind. This explanation may also apply to the SMC 
case. They further suggest that the apparent lack of slow ro- 
tating evolved stars is a spurious result, caused by erroneously 
assigning turbulent broadening - which is more pronounced for 
evolved stars (Rvans et al. 2 0021) - to rotational broadening. As 



a result, the derived v r sin i are overestimates, therefore some 
v t sin i larger than 50 km s _1 (causing the steep gradient of the 
dashed curve in the left panel at about 70-80 km s~ 1 ) in reality 
reflect projected rotational velocities below 50 km s . For this 
to be the correct explanation, the required turbulent velocities 
should be large (order 40 - 70 km s _1 ). Whether such an expla- 
nation is valid for our sample is doubtful, as we have found no 
indication for the existence of significant macroturbulent ve- 
locities in the objects with u T sin i < 100 km s . 

In the right panel of Fig.|5]the 21 unevolved SMC objects 
are compared to 66 unevolved Galactic stars as measured by 
Pennv ( 1996). Again using the Kolmogorov-Smirnov test we 
determined that the probability that the two samples have the 
same underlying distribution is 13%. Therefore, we tentatively 
assume that the SMC distribution of unevolved stars is signif- 
icantly different from the Galactic distribution. However, this 
should be verified using a larger sample of SMC objects. 

A marked difference between the two curves is the 
behaviour between the intersection points at ~90 and 
-190 kms -1 . The fact that these curves intersect at these two 
points reflects that for both the SMC and Galaxy the same 
fraction of stars show projected rotational velocities in this 
range. However, the fact that in between 90 and 190 kms 1 
the Galactic curve lies above the SMC one implies that the 
Galactic stars show preferably lower v T sin i (i.e. closer to 
90 kms -1 ) than in the SMC. This behaviour is consistent with 
the SMC stars suffering less from spin down due to mass loss, 
as the SMC stars are expected to have weaker winds. The 
behaviour outside of the above velocity range cannot be un- 
derstood within the context of spin down through winds. At 
v r sin i < 90 km s _1 the SMC shows preferably larger projected 
rotational velocities; at u T sin i > 190 km s _1 both galaxies show 
the same distribution. If the Galactic objects indeed suffer from 
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Fig. 10. Left panel: cdfs of observed (solid line) and theoretical initial u r sin i distributions. Dotted and dashed lines, respectively, 
correspond to a constant v r distribution, i.e. block function, and a single value v r distribution. Right panel: filled areas indicate 
the one <x probability distribution for theoretical u r sin i cdfs for the limited number of 21 objects calculated from a constant y, 
distribution with u m i n = Okms -1 and v lmx of 200, 352 or 600 kms~'. Shown with a dotted, dashed and dashed-dotted curve are, 
respectively, the best fitting block function, Gaussian and Maxwellian distribution cdf. See text and Fig.EJfor further comments 
and explanation. 



a stronger spin down, the latter behaviour could in principle 
imply that star formation in a relatively metal rich environment 
results in larger initial rotational velocities. However, note that 
the SMC cdf for v, sin/ > 200kms~' only contains four ob- 
jects, which makes this last statement very uncertain. 

The scenario in which SMC stars suffer less from spin 
down is consistent with the recent analysis performed by 
iDufton et ail J2006h for the young Galactic cluster NGC 661 1 . 
They find that the v t distribution of the O-type stars can be char- 
acterised by a Gaussian with a mean of ~ 125 kms _I . In the 
next section we will show that the underlying v T distribution 
of our SMC stars can be fitted by a Gaussian with a mean of 
~160kms _1 . Consequently, as the age of NGC 6611 (~2 Myr) 
is comparable to the age of NGC 346, it appears that the weaker 
winds of the SMC objects result in less spin down in the first 
few million years of evolution. 

It would be interesting to compare the evolved SMC and 
Galactic cdfs as the effect metallicity has on stellar rotation 
would likely be much more pronounced. However, for our sam- 
ple thi s is not possi ble. The reason is that in the Galactic sam- 
ple of P ennvl i 19961) the ratio of luminosity class I-II to class 
III objects is a factor two larger compared to our SMC sam- 
ple. Consequently, when making this comparison we could be 
confusing evolutionary effects with metallicity effects. 

5. Parent rotational velocity distribution 

The effect of rotation on the evolution of massive stars has been 
■ I i i I I 

studied by e.g. Heser & Lanaer (2000), Mevnet & Maeder 

(2000) and Maeder & Meynet (2001). This has shown that ro- 
tation may cause extensive mixing, changing the size of the 
reservoir of nuclear fuel availab le for evolution, a nd thus the 
lifetimes and tracks (e.g. lLanger & Maederll 19951) . The effect 



of metal content on the time evolution of the surface equa- 
torial velocity as a function of metal co ntent i s nicely illus- 
trated in Figure 10 of Mevnet & Maeder (2000) and Figure 3 
in Maede r & Mevnet! d200ll) . The first figure shows this evo- 
lution for different stellar masses in a Galactic environment, 
using an initial surface equatorial velocity of 300 kms~'. The 
second figure is identical, but now for a SMC environment. 
Ignoring the first ~10 5 yrs - which show a rapid decrease of 
Veq from ~300 to ~250 kms -1 , reflecting the initial conver- 
gence of the rotation law - the main sequence phase shows 
a monotonic decline of the equatorial rotational velocity as a 
result of an increasing radius and mass loss. Using the 20 Mq 
track as an example, u eq has reduced from -250 to ~120 km s _1 
at the end of the main sequence phase in the case of Galactic 
stars, but to only 200 kms _I for SMC stars. The main fraction 
of the SMC decline occurs near the end of the main sequence 
phase; halfway the main sequence, u eq ~ 230 km s -1 . The rea- 
son for the modest decline of the SMC star is the fact that its 
stellar wind is weaker than that of its Galactic counterparts. 
Note that for an initially 60 Mq SMC star, which has a stronger 
wind than the 20 Mq star, wind effects do play an important 
role. Most stars in our unevolved SMC sample, however, have 
masses of about ~15-30 Mq, therefore are representative for 
the discussed case. Given that the age of the NGC 346 cluster 
is about 1-3 Myr (see Sect. [7}, we may conclude that - based 
on the evolution models - the observed cdf of the unevolved 
SMC objects should lie close to the initial cdf. This allows to 
address the interesting question: what is the initial u r sin i dis- 
tribution? 
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5.1. initial v T distributions 

We compare the cdf of the unevolved SMC objects with the 
theoretical predictions. To calculate these theoretical cdfs we 
constructed simple models of the underlying unprojected rota- 
tional velocity distribution. For each distribution function we 
synthesised a theoretical distribution using a large number of 
objects (N = 10 4 ), while assuming randomly oriented rotation 
axes. The four adopted models are: 

1. a delta function, i.e. one single possible value of u r for all 
objects; 

2. a block function, i.e. a constant distribution of objects be- 
tween a minimum velocity and a maximum velocity 

^max? 

3. a Gaussian distribution, with a mean velocity v® and a stan- 
dard deviation cr; 

4. a Maxwellian distribution, specified by the most likely ve- 
locity v™; 

The normalised distributions of (best fits of the) latter three 
models are shown in the right panel of Fig. ^3 Results for 
the delta and block function are presented in the left panel 
of Fig. [TO] The top, middle, and bottom dashed curve corre- 
spond to a population with a single rotational velocity u r = 200, 
400, and 600 km s" 1 , respectively. From a comparison with the 
SMC distribution, we find that models assuming a single val- 
ued underlying u r distribution fail to reproduce the observed 
cdf. Consequently, it is very unlikely that the velocity distri- 
bution of massive stars in this cluster can be characterised by 
a single rotational velocity. This immediately implies that the 
initial v T distribution of this young cluster is not such that, for 
instance, massive stars are all born rotating at critical veloc- 
ity. Such a supposition would not be completely unreasonable 
as the initial angular momentum of Jeans-unstable molecular 
cloud fragments is so large that - if one assumes it to remain 
conserved during collapse leading to the formation of a single 
object - all stars would have to rotate super-critical. Such rapid 
rotation must lead to the ejection of material and angular mo- 
mentum, at least until the star is rotating at break-up velocity. 
As pointed out, the observed v r sin di stribution does not sup- 
port such a scenario (also see lHerrero & Naiarrol2005h . 

The dotted curves in the left panel of Fig. 1101 correspond 
to the theoretical cdfs calculated for the second model, i.e. the 
block function. The best fit to the observations requires u m [ n = 
kms 1 and o max = 352 kms -1 . For comparison, we also plot 
the result for u max = 200 (top curve) and 600 km s _1 (bottom 
curve). We see that the overall shape of this model is in better 
agreement with the observed cdf and that the best fit model 
gives a good fit to the SMC stars. 

Given the fact that the SMC cdf is constructed from a lim- 
ited number of objects, we should also try to account for the 
effect that a small sample size has on the distribution in the the- 
oretical cdfs. To do this we use the following approach. Instead 
of using a large number of simulated objects to calculate the 
theoretical cdf, we use a number equal to the amount of ob- 
served objects. Using different sets of random inclination an- 
gles we then calculate an ensemble of theoretical cdfs. The 
resulting distribution of these cdfs in the ensemble then de- 
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Fig. 11. Best fitting model v, distributions for the v t smi dis- 
tribution of the unevolved SMC objects. Shown with a dotted, 
dashed and dashed-dotted line are, respectively, the constant, 
the Gaussian and the Maxwellian distribution. 



scribes the effect of statistical fluctuations due to a limited sam- 
ple size. In the right panel of Fig.llOlthe results of this approach 
are shown for the theoretical cdfs with underlying constant v T 
distributions. The filled areas in this panel correspond to the 
ranges in the diagram containing one cr, i.e. 68 percent, of the 
theoretical cdfs. In other words the surfaces correspond to the 
area in which one can expect a theoretical cdf to be located 
within a one cr probability. The top, middle and bottom areas, 
again, correspond to distributions with y max equal to, respec- 
tively, 200, 352, and 600 km s _1 . 

So, may we conclude that the underlying rotational velocity 
distribution is indeed a block function? This would be prema- 
ture. To illustrate this we show two additional models in the 
right panel of Fig. EH The dashed curve is the best fit for a 
Gaussian u r distribution, with a mean velocity - 158 kms -1 
and cr = 97 kms" 1 . The dashed-dotted curve is the best fit for 
a Maxwellian distribution, characterised by a most likely ve- 
locity v^f = 148 kms -1 . For reference, the dotted line gives 
the best fit for the constant distribution. The Gaussian and 
Maxwellian models have an uncertainty which is similar to 
that of the block function (the dotted line, with the middle grey 
area representing the uncertainty). Given the modest difference 
between the block and the Gaussian curves - and in view of 
their error bars - it is clear that we cannot distinguish between 
these two models. The Maxwellian function is also very simi- 
lar, however, here the agreement with the observations seems a 
bit off, in particular at low velocities. 

Why are the block and Gaussian distribution so indistin- 
guishable? The reason must be that the underlying intrinsic 
distributions are essentially similar. This is shown in Fig. ^2 
Though the models obviously show differences, they are both 
characterised by a mean velocity of about 150-180 kms -1 and 
a comparable effective half width of roughly 100-150 km s _1 . 
Therefore, this is currently the most meaningful and robust 
specification we can assign to the underlying rotational veloc- 
ity distribution. The similarity between the models in the cdf 
plot shows that even with a significantly larger set of observed 
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rotational velocities it will be difficult to better define the exact 
shape of the parent population of v r . The Maxwellian distri- 
bution is a bit worse, essentially because it is lacking in stars 
with low (< 50 kms~' say) rotational velocities. We conclude 
that at birth the massive star population of NGC 346 must have 
included ~5— 15 percent slow rotators. 

6. The weak wind problem 

The majority of our programme stars have luminosities less 
than L* ~ 10 5 25 Lq. The winds of these SMC O- and B-type 
stars are currently at the focus of attention as they do not ap- 
pear to agree with theoretical predictions JBouret et alJ feoCB: 
iMartins et all2 004). Starting from this luminosity the observed 
rates rapidly fall below the predicted ones, leading to a discrep- 
ancy of about two orders of magnitude at L* ~ 10 5 Lq. The 
reason for this discrepancy is not understood so far. The cause 
may be ill treated or missing physics in the M predictions, e.g. 
because of a break down of the adopted Sobolev approxima- 
tion for low density winds JOwocki & Pulsll 1 9991) o r possibly 
the neglect of ion-decoupling ( Kitic kTetalTl2003l) . Whether 
the last hypothesis is indeed a viable option seems doubtful. 
Computations of Krtic kaet alJ appear to indicate that this ef- 
fect only starts to play a role a t a metallicity Z < 1/30 Zq. 
Moreover. lMartins et alJj200 5b) also report a weak wind prob- 
lem for Ga lactic stars m aking ion-decoupling even more un- 
likely. The IMartins et alJ result also suggests that metallicity 
effects can be excluded. An alternative explanation may be that 
the dwarfs showing the weak wind problem represent an ear- 
lier evolutionary state, i.e. a state in which the wind is not yet 
"fully developed". The problem with this hypothesis is that the 
characteristic timescale for the wind to develop - of order of 
the dynamical timescale of the wind - is very short compared 
to evolutionary timescales, and that the wind properties are de- 
termined by global and atmospheric properties only. 

The cause of the discrepancy may also be connected to a 
defect in the spectroscopic derivation of the mass loss rate. 
The winds of these SMC dwarf O-type stars are quite weak 
(~10~ 8 - lO~ 6 M0yr~'), and reach the limit s of sens itivity of 
Ha as a mass loss diagnostic. Therefore, the Boure t*et"afl and 
IMartins et alJ study use unsaturated ultraviolet resonance lines 
of carbon, nitrogen, and oxygen, which are much more sensi- 
tive mass loss indicators. However, these lines typically belong 
to trace ions, for which the ionisation is extremely sensitive to 
the local (shock generated) X-ray radiation field. Though sig- 
nificant progress in our understanding of the processes lead- 
ing to non-thermal X-rays in stellar winds has been gained in 
the last years (e.g. Feldmeier et al. ll997lJ Pauldracfi et alj200lt 
iRramer et all2 003). we cannot yet exclude problems in the UV 
based M determination as the cause for the weak wind discrep- 
ancy. 

Our genetic algorithm based fitting method allows to de- 
termine the mass loss at the "upper end" of the weak wind 
regime. Figure[Q]shows a comparison of six O stars for which 
the UV and, save for NGC 346-026 and NGC 346-028, optical 
spectra have been used to determine M JCrowfher et alJ l2002: 
iBouret et all2003l) and which we have analysed here using op- 
tical spectra only. The four stars with the highest mass loss all 
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Fig. 12. Comparison of mass loss rates determined in this study 
using the optical spectrum with values determined from the UV 
spectrum. For the three objects with the lowest mass loss rates 
the values determined from the optical spectra correspond to 
upper limits. 

have a luminosity L > 10 5 3 L©, i.e. they are in a regime in 
which the theory agrees with observations. In these cases the 
UV+optical studies all include Ha. The M results reported in 
these studies compare very well with our findings. 

The two remaining stars, NGC 346-026 and NGC 346-028 
have been analysed by Bou ret et alJ (12003). Their luminosities 
are 10 4 93 and 10 5 10 Lq, respectively. This is in the regime 
wh ere observed and predicted mass loss disagree. Though 
fhe lBouret et alJ analysis includes part of the optical spectrum 
(from 3910 to 5170 A), it does not include Ha. For NGC 346- 
026 they derive M = 3.2 x 10" 9 M Q yr _1 ; we find a -30 times 
higher rate. For NGC 346-028 they obtain M = lO^Moyr -1 , 
where we find a ~500 times larger value. Note, however, that 
for both stars the downward error bars on our fit values ex- 
tend to the lower limit of the mass loss regime in which the 
automated method was allowed to search for a solution. This 
implies that our results are de facto upper limits. Therefore, we 
cannot (yet) conclude whether or not the UV and Ha diagnostic 
agree or disagree in weak wind situations. 

Note that though we interpret these M as upper limits, the 
automated fitting method did return a best fit value. Figure [8] 
shows the strength of the stellar winds of our entire sample 
in terms of the modified wind momentum (see Sect. 14.61 for a 
discussion). The inverted triangles represent upper limits. All 
these upper limits are well above (up to two orders of mag- 
nitude) the mass loss rates derived on the basis of UV profile 
fitting; not a single one "by chance" gives a value comparable 
to what is found with the UV method. Why is this so? We can 
identify two potential reasons. First, the automated fitting pro- 
cedure scans the mass loss dimension in a linear way. For weak 
wind cases typical lower and upper bounds of this dimension 
are 10~ 9 and 10~ 7 MQyr~'. In most cases it is reasonable to 
assume that only above 10~ 7 Moyr~' the profile of Ha is sen- 
sitive to mass loss. In that case it is, e.g. nine times as likely that 
the automated method more-or-less by chance settles for a mass 
loss above 10~ 8 MQyr~\ and not for the much lower values in- 
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Fig. 13. Comparison of the programme stars locat ed in NGC 346 with iso chrones for zero age, 1, 2, 3, 4, 5 and 7 Myr de- 
rived from the non rotating evolut ionary models of Charbonne l et alJ (119931) (dashed lines) and models including rotation from 
Maeder & Mevnet (2001) and lMevnet & MaedeJ J2005i) (dotted lines'). For reference the tracks of the rotating models are also 
shown (grey lines, see Fig.[6]for initial masses of these tracks), where for clarity purpose only the blueward part of the evolution 
of the most massive tracks is plotted. Shown with identical symbols as in Fig.[5]are our programme stars, which are located in 
NGC 346. Luminosity class Vz objects are denoted using grey circles. See text for a discussion. 



dicated by UV analysis. Second, small errors in the broadening 
function of Ha may perhaps be "corrected" in the automated 
fitting process by settling for a mass loss which is "as high as 
possible". Consequently, we conclude that we cannot assign 
any importance to the best fit value returned by the automated 
method when the downward errors extend to the lower bound 
of M allowed in the fitting process. 



7. On the evolutionary status of NGC 346 

Undoubtedly the large number of early-type stars in NGC 346 
dMassev et all 1989) is indicative for a young age. Further evi- 
dence for the youthful nature of this cluster may be given by the 
presence of several so-called Vz stars, which are hypothesised 
to represent the e arliest stages of main sequence evolution o f 
massive stars (e.g. Walborn & Parker 1992; Parker et al. 
In this section we will determine an age estimate for NGC 346 
and try to quantify the uncertainty rotation introduces in this 
estimate. We will also discuss how well the spectroscopic Vz 
designation correlates with ZAMS evolution. 



7.1. The age of NGC 346 

Rotation affects both the lifetimes and the tracks of stellar 
evolution models. Consequently, before we can estimate the 
age of NGC 346 we need to assess the systematic uncertainty 
rotation will introduce in the age estimate. We do this by 
considering the difference s between isochrones de rived from 
the non rotating models of Charbonnel et al. ( 1993) and those 
based on the rotating mod els of lMaeder & MevneJ 12001 ) and 
iMevnet & MaedeJ d2005h adopting an initial v T of 300 km s -1 . 
In Fig.^]the two sets are shown as, respectively, dashed and 
dotted curves for zero age, 1, 2, 3, 4, 5 and 7 Myr. For reference 
the evolutionary tracks of the models accounting for rotation 
are also shown in this figure using grey lines. 

Age determinations of individual stars, based on the evo- 
lutionary tracks, are given in Fig.[2] The error bars on these 
life times account for the uncertainties in temperature and lu- 
minosity. Open symbols (see the figure caption for details) re- 
fer to non-rotating tracks; grey symbols for those accounting 
for rotation. Note that they differ only very modestly for rel- 
atively unevolved objects. For instance, stars near the main 
sequence are judged to be ~10 5 years younger if rotation is 
taken into account. We estimate that for main sequence stars 
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Fig. 14. Age determination of the programme stars in NGC 346 
using the non-rotating SMC models of ICharbonnel e t al. 
( 199 3b (open symbols' ) and those including rotation 
1 Maeder & Mevnetl l200l( iMevnet & Maederl 120051) (closed 
symbols). Circles, triangles, and squares denote dwarfs, giants, 
and supergiants, respectively. The horizontal axis gives the ID 
number of the star (see e.g. Tab. [3}; the vertical axis the age in 
Myr. The lower limits (upwards pointing arrows) provide age 
estimates from tracks of chemically homogeneous evolution 
(Yoon et al. in preparation); see text for a discussion. 



having luminosities larger than logL+/Lo « 4.5 this system- 
atic difference may increase up to approximately 0.5 Myr. For 
the more evolved phases, i.e. approximately giant phase and 
later, this systematic discrepancy switches sign and rapidly 
becomes larger. Considering the supergiant NGC 346-001 at 
log L+/Lq = 6.0 we estimate its age to be ~1 Myr more if it 
is initially rotating with u T ~ 300km s _1 compared to the non- 
rotating case. 

Ignoring the supergiant NGC 346-012 at log off » 4.4 
for w hich cluster membership is very uncertain JEvans et alJ 
2006), we find that the objects populate the region between the 
ZAMS and the 7 Myr isochrone, with typical error bars (in both 
directions) of 2-5 Myr. 

So, what could explain the apparent age scatter? First, some 
objects in our sample might actually not belong to the cluster. 
In addition to NGC 346- 012 this is, for instance, also suggested 
bv lWalborn et"aTI (EoOO) for the B0 dwarf NGC 346-026. Based 
on its spatial location, discrepant radial velocity and stellar pa- 
rameters these authors conclude that this object is not a coeval 
member of the cluster. Its large age of about 7 Myr seems to 
support this conclusion. Because of cluster membership issues 
we have placed both NGC 346-012 and NGC 346-026 to the 
right in Fig. [H] and have separated them from the other stars 
using a vertical dashed line. For the remainder of the discussion 
we treat them as not belonging to the cluster. Second, the Oe 
star NGC 346-018 is poorly fitted because of contamination by 
circumstellar material. This implies that the parameters derived 
for this star should be taken with considerable care, and there- 
fore also the apparently large age of about 7 Myr that we derive 
for this object. Ignoring the above three discussed objects the 
oldest investigated star in NGC 346 is ~ 5 Myr (NGC 346-010). 



What remains is an age spread of stars from zero to the 
above motivated 5 Myr. Can such an age spread be explained? 
We first checked whether there is a relation between the age of 
the stars and their distance to the cluster centre. Though we do 
see that the oldest stars (NGC 346-001 through NGC 346-025) 
are all at some distance from the centre this also holds for some 
that are assigned very young ages. Also two relatively old stars 
(NGC 346-007 and NGC 346-022) are located relatively close 
to the cluster core. Therefore, we can not claim such a relation. 

Next we consider the possibility that the adopted tracks are 
not appropriate for some stars. As we have explained, when 
it comes to age determination there is no fundamental differ- 
ence between tracks without rotation and those that do account 
for it starting the evolution from an initial v T = 300 kms -1 . 
However, as already discussed in Sect. 14.41 if u r exceeds some 
critical value the evolution will show a bifurcation behaviour, 
leading to tracks that evolve bluewards (from the ZAMS) and 
upwards. This is the result of such efficient mixing that the 
HRD behaviour proceeds as for chemically homogeneous evo- 
lution. Adopti ng such tr acks from Yoon et al. (in prep.) (also 
see lYoon & Lange3l2005l) it turned out to be difficult to de- 
rive ages as for near critical rotation the effective temperature 
changes so rapidly (as a function of v r ) that one can in practise 
only rely on luminosity as an age discriminator. Given the error 
bars on L and M s this did not lead to meaningful constraints for 
most of the objects. Moreover, in a few cases the mass discrep- 
ancy turned out to be so large that the luminosity range implied 
by M s (plus errors) did not coincide with the observed L+. This 
also prevents one from deriving an age estimate. 

A different approach we tried is based on the fact that for 
chemically homogeneous evolution the surface helium abun- 
dance is essentially a measure of stellar age, for given initial 
mass. Taking into account the errors on Yn e and M s , this may 
yield lower limits to the age of the stars. One derives lower lim- 
its because if evolution is not fully chemically homogeneous 
the measured helium abundance is a lower limit for the mean 
stellar Fee- The method will only provide useful results for 
those stars that show enhanced helium even when correcting 
for the negative error in Fee- This is the case for four objects 
(NGC 346-001, NGC 346-028, NGC 346-31, and NGC 346- 
050). The derived lower age limits are shown in Fig.ll4l(arrow 
symbols). The lower limit of NGC 346-001 (the most massive 
star) is consistent with the age as derived using the standard 
method. Interestingly, the lower limits of the ages of the three 
remaining stars are all significantly above what we derived us- 
ing modestly rotating or non-rotating models. We have no indi- 
cation that our Fee determinations are flawed or that these stars 
had an initial helium fraction significantly higher than 0.09. 
Consequently, it is likely that these stars evolve along chem- 
ically homogeneous tracks and have ages comparable to the 
oldest stars in the cluster. 

In view of the above discussion we have to conclude that 
we can not derive a unique age for NGC 346. The distribution 
of ages seems to indicate two preferred values: 1-3 Myr and 
3-5 Myr. We add that some stars appear to reside on or even 
to the left of the ZAMS. Though, it could be that these lat- 
ter stars evolve along tracks for rapidly rotating stars (follow- 
ing tracks for chemically homogeneously evolving stars). If so, 
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we could not constrain their ages from such models. The ages 
estimated from modestly rotating or non-rotating evolutionary 
models of three helium enriched dwarf objects are so low (less 
than 2 Myr) that it is not possible to bring them into agreement 
with the derived Yu e values. Chemically homogeneous evolu- 
tionary models instead predict an age comparible with the old- 
est cluster population. Given/assuming that all stars belong to 
the same compact cluster it seems very unlikely that the range 
of ages implies that NGC 346 experienced a prolonged burst of 
star formation, or a series of short bursts resulting in multiple 
coeval populations. 

Th e age of NGC 34 6 was a ls o estimated bylKudritzki et alJ 
ll 1983b. IMassev et alJ i 1 9891). IWalborn et alJ l l2000l) and 
Bou ret et alJ J2003I) ("also see lMassev et al.ll2.005h . In the first, 
second and fourth studies ages of 2.6 Myr, 2-4 Myr and 3 Myr 
were found, w hich in principle are consistent with our findings. 
IWalborn et alJ estimate an age of only 1 Myr. This agrees with 
our lower age limit. The reason for their young age estimate is 
that their four star NGC 346 sample predominantly contained 
the stars in the cluster c losest to the ZAMS. Combined with 
the T e ff calibration from lVacca et alJ \ 19961) adopted by these 
authors, which resulted in relatively high estimate s for this pa- 
ramete r, for instance, compared to the analysis of B ouret et alJ 
J2003h . explains their relatively young age estimate. 

7.2. ZAMS stars 

The defining characteristic of the luminosity class designation 
Vz is a very intense He n A4686 absorption line that is substan- 
tially stronger than any of the other helium absorption lines. 
The distinct behaviour of this line is believed to be linked to a 
luminosity effect connected to an early evolutionary phase. As 
Vz stars are hypothesised to lie relatively close to the ZAMS 
they are anticipated to be less luminous compared to normal 
dwarf stars. As a result of this, the relatively strong /14686 line 
can be explained as a lack of partial filling in of the line profile 
due to (wind) emission which normal dwarf stars do experience 
JParker et alii 992HWalborn & Bladesll 9971). 

In the HR-diagram shown in Fig. ^]we have highlighted 
the location of the three NGC 346 Vz stars in our sample 
using grey circles. From these three objects, NGC 346-028, 
NGC 346-031 and NGC 346-051, we find that the first and last 
are indeed found to lie on or left of the ZAMS. NGC 346-03 1 in 
contrast occupies a location between the one and two million 
year isochrones. We suspect that that the strong He n /14686 
line in the spectrum of this object is in fact related to its pecu- 
liar chemical composition (see Sect. l4.4> . Moreover, the spatial 
separation of NGC 346-03 1 of ~5 arcminutes from the cluster 
centre, makes it unlikely that this object is formed recently as 
part of the cluster. Consequently, care should be taken with the 
correct interpretation of Vz stars as near-ZAMS objects, based 
on optical classification criteria only. 

In addition to the two Vz stars Fig.[0]reveals the presence 
of three more stars located on or left of the ZAMS. Even though 
the spectra of these objects, NGC 346-046, NGC 346-107 and 
NGC 346-112, do not exhibit the Vz characteristics, their lo- 
cation in the HR-diagram clearly suggest an early evolutionary 



phase. Why is this so? We suspect that the reason for the lack 
of a strong He n /I4686 line in the spectra of these objects is re- 
lated to their relatively low effective temperatures. Compared 
to the neutral helium lines the relative strength of the He n lines 
decreases for decreasing T e g, Consequently, even though its 
line profile is not filled in it becomes more unlikely for the 
He ii /14686 to be stronger than all the He i lines for relatively 
low effective temperatures. In Fig.^]we see that these normal 
dwarf stars are indeed cooler than the Vz stars at the ZAMS. 
Interestingly, it turns out that in the spectrum of the hottest 
of these three objects, NGC 346-046, the He n A46S6 is the 
strongest He n line, which has a strength that is approximately 
equal to the strength of the strongest He i line. We therefore ar- 
gue that the effective temperature of this object, i.e. ~40 kK, 
corresponds to the lower limit of T e g for SMC metallicity for 
which a Vz luminosity class designation is at all possible. 

8. Summary and conclusions 

Using the automated fitting method developed by 
iMokiem et all J2QP5), which combines the stellar atmo- 
sphere code fastwind with the genetic algorithm based 
optimisation routine pikaia, we have performed a quantitative 
analysis of 31 O- and early B-type stars located in the SMC. 
This sample was mostly drawn from the targets observed 
in the SMC clusters NGC 346 and NGC 330 as part of th e 
VLT-FLAMES survey of massive stars (see lEvans et all20 05). 
and is the largest sample of early-type SMC stars analysed so 
far. Even though many of the observed spectra show nebular 
contamination and have signal-to-noise ratios as low as 50, 
the fitting method did not encounter convergence problems. 
Instead we find that the quality of the data is naturally reflected 
in the errors that are estimated for the fit parameters using the 
width of the global optimum. 

Concerning the stellar properties of the objects in our sam- 
ple we can draw the following conclusions: 

i) The average effective temperature for a given spectral type 
for SMC dwarf stars is hotter by typi cally 3 t o 4 kK than 
for their Galactic counterparts (Marti ns et alJ l2005a). We 
attribute this hotter temperature sc ale to a reduced line blan- 
keting in the SMC. In contrast to IMassev et al 1 (120051) we 
find that this blanketing effect on T e s extends to spectral 
types as late as B0. 

ii) The helium enrichment found in our stars does not appear 
to be consistent with predictions of the evolution of mas- 
sive stars ignoring rotation. Comparison with evolution- 
ary tracks that do account for rotation shows a qualitative 
agreement. In the HR-diagram the region in which a he- 
lium enrichment is predicted roughly coincides with the 
location of evolved objects showing evidence of Tne en- 
richment. The actual degree of enrichment is in many cases 
still strongly under predicted by the most recent evolution- 
ary models. 

Hi) From our three programme stars with a Vz luminosity class 
designation, we find that only two are true ZAMS objects. 
Moreover, three additional objects, which are not classified 
as Vz stars, are found to lie on or left of the ZAMS. We 
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argue that this is related to a temperature effect inhibiting 
relatively cool stars located on the ZAMS to exhibit the 
spectral features characteristic for the Vz luminosity class. 

iv) We find good agreement between our spectroscopically de- 
rived masses and those derived from evolutionary calcula- 
tions. A mild mass discrepancy is still found for stars with 
M s < 20 Mq. This discrepancy correlates with the surface 
helium abundance and is qualitatively consistent with the 
predictions of chemically homogeneous evolution. 

The stars in the young cluster NGC 346 allow to study the 
role of mass loss and angular momentum loss during the early 
evolution of massive stars in a low metallicity environment. 
Regarding these issues we can conclude that: 

v) The observed distribution of ^ sin / for evolved stars (lu- 
minosity classes I-II) contains relatively less fast rotators 
and slow rotators compared to the distribution for un- 
evolved stars (luminosity classes IV-V). These findings are 
in agreement with analyses of Galactic samples and can 
be interpreted as a spin down due to an increased radius 
and angular momentum loss through the stellar wind, and 
excess turbulent bro adening among the evolved stars (cf. 
iHowarth et alJ 19971) . We also find that compared to the ve- 
locity distribution of unevolved Galactic objects the distri- 
bution for unevolved SMC objects shows a relative excess 
of stars rotating with projected velocities 120 > i> r sku <, 
190 km s -1 . This excess can be interpreted as evidence of 
a reduction of dissipation of angular momentum in weaker 
stellar winds resulting from the low Z environment of the 
SMC. 

vi) Using a simple modelling approach we have constrained 
the underlying v, distribution of the unevolved SMC stars. 
The exact shape of this underlying distribution - for which 
we have tried a block, Gaussian and Maxwellian shape - is 
found to be degenerate, however, it can be characterised by 
a mean velocity of about 150 - 180 km s _1 and an effective 
half width of roughly 100 - 150 kms -1 . The Maxwellian 
distribution appears a somewhat poorer representation of 
the v, distribution, indicating that 5-15 percent of the stars 
must be slow rotators (u r < 50 km s _1 ). 

vii) It is not possible to derive a unique age for NGC 346 us- 
ing moderately rotating or non-rotating evolutionary mod- 
els. Though the bulk of the stars indicate a lifetime of 1- 
3 Myr, part of the investigated stars seem somewhat older 
(3-5 Myr). Clearly the ages depend on the adopted evolu- 
tionary tracks. However, we point out that predictions for 
non-rotating stars and those accounting for the rotational 
evolution of stars with an initial u r = 300 km s _1 essentially 
give the same result. Three main sequence stars, all with 
an age determination of less than 2 Myr, show significant 
helium enhancement. Assuming that these stars have the 
same initial Yu e this implies a lower limit to their age of 2 
to 3.5 Myr. 

viii) Based on the predictions of the time evoluti on of the equa- 
torial velocity for low metallicity ( Maeder & Mevnet 2001 ) 
the above age considerations imply that the underlying u r 
distribution derived for our sample (conclusion vi) is repre- 
sentative for the initial rotational velocity distribution. 



ix) Mass loss rates determined for the SMC objects are found 
to be systematically lower than for Galactic objects. For 
log L+/Lq > 5.4 we find that the wind strengths are 
in excellent agreement w i th the theoretical predictions of 
Wink, de Koter. & Lamers\ \200l\) . The empirical modified 
wind momentum luminosity relation constructed for this 
luminosity regime agrees to within the error bars with the 
predicted WLR. 

x) Down to log L+/Lq w 5.3 the mass loss rates determined 
from the optical spectra are in agreement with UV rates. 
For luminosities lower than this value the winds become 
too weak for an accurate determination of M from the opti- 
cal spectrum only. Consequently, we cannot (yet) conclude 
whether or not the UV and Ha diagnostics agree or disagree 
in weak wind situations. 
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Appendix A: Fits and comments on individual 
objects 

The observed spectra shown in this section were corrected for 
radial velocities. If not noted differently the lines that were fit- 
ted are the hydrogen Balmer lines Ha, Hy and H6; the Hei 
singlet line at 4387 A; the Hei triplet lines at 4026, 4471 and 
4713 A; and the He n lines at 4200, 4541 and 4686 A. In general 
we find that all lines could be fitted quite accurately including 
the He i AAAll for which we only encountered the "generalised 
dilution" in one case. This is in accordance with the fact that 
the majority of the sample consists of dwarf stars and the sin- 
gle star where this problem is present, AzV 372, is a late type 
supergiant. 

NGC 346-001 The fit of the spectrum of this 07 supergiant 
is presented in Fig lA.ll All lines except Hei /14387 are repro- 
duced well. Note that the nebular contamination on top of the 
Ha line core was disregarded in the fit. The inaccurate repro- 
duction of Hei /14387 is the result of t he low w eight given to 
this fine in our weighting scheme (see IPaper J) . For this low 
weight the quantitative fit quality of Hei /14387 is comparable 
to that of Hei AAQ26. Therefore, the discrepancy is not signifi- 
cant and merely reflects the uncertainty with which we believe 
that the stellar atmosphere code can reproduce this line. 

Th is object was previously studied by JCrowther et al. 
( 2002) using the alternative wind code cmfgen dHillier & M iller 
1998). We find good agreement between the photospheric pa- 
rameters. Note that the helium abundance Tee = 0.2 adopted by 
ICrowther et al.l is in good agreement with the F He = 0.24 self 
consistently determined here. The largest difference is found 
for the effective temperature (2 kK). This difference can be 
explained by the higher surface gravity obtained from our fit, 
which requires a higher to obtain the same helium ioni- 
sation equilibrium. This occurs in many of the fits discussed 
here, and we attribute this to a better multidimensional param- 
eter optimisation by the automated method. This assures that 
the global optimum in parameter space is found avoiding log g 
and T e s val ues whic h could correspond to a local optimum (see 
Sect. 5.2 of lPaper J) . 

The mass loss rate obtained from our fit of the optical spec- 
trum is within the error bars in agreement with the value that 
ICrowther et al J derived from the UV spectrum. With respect to 
the acceleration of the wind we find a large differ ence. From the 
optical spectrum we derive /? = 1.15, whereas fCrowther et alJ 
find p — 1.65 based on the UV resonance lines. They do note 
that from their modelling of Ha they obtain a smaller value 
for this parameter, indicating a possible systematic discrepancy 
between wind accelerat ions dete rmined from optical and UV 
spectra (also see lPuls et all 1 996). 

NGC 346-007 With exception of the Hei line at 4471 A, 
which suffers from very strong nebular contamination, a good 
fit was obtained for this object. Note that the fit shown in 
Fig. lA.ll does not include the Hei 4713 line. The nebular con- 
tamination in this line proved to be too severe to obtain a reli- 
able fit. 



Within the uncertainties the photospheric pa rameters T e jf 
and l ogg are in good agreement with those of iBouret et alJ 
(2003) who used cmfgen for their analysis. These authors find 
T e ff = 41 .5 kK and log g = 4.0. The automated method was not 
able to self consistently derive /?. Consequently, a fixed value 
of /3 = 0.8 was used. With this fixed value a mass loss rate of 
2.3 x 10~ 7 Moyr -1 is f ound, which i s also in good agreement 
with the rate derived by Bo uret et alJ from a fit to the wind sen- 
sitive UV lines. 

The effective temperature of T e g = 42.8 kK obtained here 
is in poor agreement with the recent analysis of iHeap et alJ 
(2006). Using the plane-parallel atmosphere model tlusty 
these authors derive a relatively low r e jf of 40 kK. The rea- 
son for this discrepancy is unknown, but can probably be at- 
tributed to a different treatment of line blanketing in tlusty 
and fastw ind, or the coarsely sampled grid of models used by 
IHeap et alJ 

NGC 346-010 is a fast rotator. To obtain the fit displayed 
in Fig. lA.ll a projected rotational velocity of 313 kms 1 was 
needed. Even though all lines seem to be fitted correctly we 
note a slight over prediction of the strength of the wings of 
the wind sensitive He n /14686 and Ha lines. The shape of the 
former line could indicate the existence of denser equatorial 
material or possible the presence of a disk. 

NGC 346-012 As can be seen in Fig. IA. li the low projected 
rotational velocity of this B-supergiant allowed for a nearly 
perfect fit. Note that given the relative scaling of the plots in 
Fig. IA. li the fit quality of the He n lines is similar to that of the 
neutral helium lines, i.e. they are fitted well. 

The good fit quality obtained resulted in relatively small 
error estimates (cf. Tab. @}. Most interestingly, the uncertain- 
ties on the wind parameters, considering the very tenuous wind 
of this object, are modest as well. This is mostly because Ha 
is not affected by nebular contamination. As a result of this 
the strength and the shape of this wind sensitive line could be 
matched very accurately with little room for error. 

NGC 346-018 The line profiles of this object suffer from 
strong contamination, making it problematic to obtain a good 
fit. We believe that this emission is not only nebular in origin, 
but, given the breadth in the hydrogen lines, is possibly also 
emitted from a circumstellar disk. The contamination in the Ha 
Balmer line proved so severe that we did not include it in the 
fitting process. Our best fitting model is presented in Fig. IA.2l 

NGC 346-022 A good quality fit was obtained for this object. 
The nebular emission in the core of the He i 4471 line did not 
hamper an accurate determination of the effective temperature, 
hence the small error given in Tab.|4] This T e ff is considerably 
larger than the effectiv e temperature one w ould expect for a 
Galactic 09 dwarf (see Mar tins et alJl2005al) . However, it is in 
line with the expectation of a reduced line blanketing effect as 
a result of a low Z environment, and is in agreement with the 
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Fig. A.l. Comparison of the observed line profiles of NGC 346-001, -007, -010 and -012 with best fitting synthetic line profiles 
obtained using the automated fitting method (grey lines). Wavelengths are given on the horizontal axis in A. The vertical axis 
gives the normalised flux. Note that this axis is scaled differently for each line. 



values obtained for the other 09 V stars in our sample (see 
Sect.l4~TV 

iHean et all §006) determined a relatively low T e g of 
32.5 kK for this object, which is a consequence of the low sur- 
face gravity of loggr = 3.75 determined by these authors. In 
contrast the automated method obtained log g - 4.2 requiring a 
higher to obtain the correct helium ionisation equilibrium. 



NGC 346-025 The parameters obtained for this object are 
very similar to the parameters used to fit the previously dis- 



cussed object. The only difference is in the projected rotational 
velocity, which is higher by more than a factor of two. 



NGC 346-026 All lines, including the weak He n lines, are fit- 
ted well for this object. In the Fig. lA.2l the final fit is presented. 
The luminosity classification IV ascribed to this object is re- 
flected in the gravity of log g = 3.76. This value is still higher 
by 0.1 6 de x compared to the value determined by B ouret etaD 
(2003]) and lHeapetaD (l2006). Accordingly, the T e $ determined 
by the automated method is also higher by 1 . 1 kK. 
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Fig. A.2. Same as Fig.|AT] however, for NGC 346-018, -022, -025 and -026. 



NGC 346-028 This object is classified as an OC6 Vz star, 
indicating that its spectrum shows signatures of a relatively 
low nitrogen abundance and evidence for an early evolutionary 
state. Evidence for the latter is provided by the large relative 
strength of the Hen i4686 line (also see Sect.l72t, In Fig. IA.3l 
the best fit is presented. Note how well the He n blends in the 
hydrogen lines are reproduced. The effective temperature and 
luminosity place this object at a position close to the ZAMS in 
the HR-diagram shown in Fig.^J corroborating its young age. 



Compa r ed to the analyses of iBouret etal] J2003I) and 
iHeap et all J2006) we find an effective temperature higher 
by almost 3 kK. This is likely related to the large helium 



abundance we found. In contrast. IBouret et alJ and lHeap et all 
adopted the Solar value for this parameter. In terms of the 
wind parameters we were only able to determine an upper 
limit of lO~ 7 M0yr~' whic h is at least consistent with the 
3 x 10~ 9 Moyr~' derived bv IBouret et alJ from the UV spec- 
trum. 



NGC 346-031 is also classified as a Vz type object. The fit to 
its spectrum is shown in Fig. IA.3l In contrast to NGC 346-028 
the location of this object in the HR-diagram does not corrob- 
orate its youthful spectral classification. Instead the isochrones 
suggest an age of 1-2 Myr. 



28 



M. R. Mokiem et al.: Mass loss and rotation of early-type stars in the SMC 




346-028 OC6 Vz 

Ha 



Hel + Hell 4026 




4325 4340 4355 
Hel 4387 



6542 6562 6582 
Hel 4471 



4018 4026 4034 
Hel 4713 




4190 4200 4210 



4534 4541 4548 



4680 4686 4692 



346-031 08 Vz 

Ho 



Hel + Hell 4026 



- 0.60 - 




4325 4340 4355 6542 6562 6582 4021 4026 4031 




0.60 - 





4386 4388 4390 4468 4471 4474 4712 4713 4714 

Hell 4200 Hell 4541 Hell 4686 

1.00 




4190 4200 4210 4534 4541 4548 




4686 4692 



346-033 08 V 

Ha 



Hel + Hell 4026 







1.00 | 






0.90 - 






0.80 - 




4325 4340 4355 6542 6562 6582 

Hel 4387 Hel 4471 



4017 4026 4035 
Hel 4713 




4381 4388 4395 4463 4471 4479 4709 4713 4717 

Hell 4200 Hell 4541 Hell 4686 

1.00 




4190 4200 4210 



4533 4541 4549 



4680 4686 4692 



346-046 07 Vn 

Ha 



Hel + Hell 4026 




0.84 - 

4325 4340 4355 6542 6562 6582 4017 4026 4035 

Hel 4387 Hel 4471 Hel 4713 




4381 4388 4395 4462 4471 4480 4705 4713 4721 



1.00 
0.95 
0.90 



5 



AM 




4190 4200 4210 



4534 4541 4548 



4680 4686 4692 



Fig. A.3. Same as Fig.|A7T] however, for NGC 346-028, -031, -033 and -046. 



NGC 346-033 The fit presented in Fig.|A3lshows that despite 
the strong nebular contamination in the neutral helium lines, a 
good fit could be obtained for this 08 dwarf. The best fit mass 
loss rate of 7.4x 10~ 7 M0yr _1 seems rather high for its luminos- 
ity. The reason for the large M is likely the anomalously large 
wind velocity (~4100 kms -1 ), which resulted from a scaling 
with the escape velocity at the stellar surface (see Sect. 14. 6> . 



NGC 346-046 The shape of the line profiles of this object 
are indicative for fast stellar rotation. This is confirmed by the 
fit shown in Fig. IA.3I for which a i> r sin/ of 340 kms -1 was 
needed, classifying 346-046 as a very fast rotator. The helium 



abundance is slightly enhanced, which could be related to the 
rapid rotation of this object. However, one has to keep in mind 
that the error estimate on this parameter is of the same order of 
the enrichment, when considering this statement. 



NGC 346-050 To fit the line profiles, displayed in Fig. IAT4I 
a u r sinz of 357 kms -1 was needed. Interestingly, similar to 
the other fast rotator 346-046 a small helium enhancement was 
found. 

NGC 346-051 is the final Vz object in our sample. Its spec- 
trum suffers from strong nebular contamination. Nevertheless 
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Fig. A.4. Same as Fig. EI] however, for NGC 346-050, -051, -066 and -077. 



a good fit could be obtained, which is presented in Fig. I A. 41 NGC 346-077 A good fit was obtained despite the strong neb- 
The location of NGC 346-051 in the HR-diagram presented in ular contamination affecting the hydrogen and neutral helium 
Fig-Elpoints to an early evolutionary phase for this object. lines. The final fit is displayed in Fig. IA.4l 



NGC 346-066 The line profiles of this object are matched well 
by the synthetic profiles. In the final fit, shown in Fig. I A. 41 
only the Hei line at 4713 A shows a slight mismatch due to 
an apparent over prediction of the width of this line. Given the 
signal-to-noise ratio of the spectrum, we do not believe this to 
be significant. 



NGC 346-090 Like 346-066 this object is of spectral type 
09 V. The fit to its spectrum is presented in Fig. IA.5l Compared 
to 346-066 a T e ff lower by ~800 K was obtained. This is the re- 
sult of the weaker He n lines. 



NGC 346-093 The fit to the spectrum of the B0 dwarf 
NGC 346-093 is shown in Fig. IA.5l This fit shows that the very 



30 



M. R. Mokiem et al.: Mass loss and rotation of early-type stars in the SMC 




346-090 09.5 V 

Ha 



Hel + Hell 4026 



I, 




4325 4340 4355 6542 6562 6582 

Hel 4387 Hel 4471 

1.00 



4017 4026 4035 
Hel 4713 





4381 4388 4395 4462 4471 4480 




4709 4713 4717 
Hell 4686 



1.00 


ilk! 
It 




1.00 


0.97 






0.95 


0.94 




f I; 


0.90 




4190 4200 4210 



4533 4541 4549 



4680 4686 4692 




346-093 B0V 

ht 



1.00 






0.90 


- 




0.80 






0.70 







4080 4100 4120 
Hel + Hell 4026 



4325 4340 4355 
Hel 4387 




401 8 4026 4034 
Hell 4200 



4382 4389 4396 
Hell 4541 



0,5 : W 




6542 6562 6582 
Hel 4471 




4463 4471 4479 
Hell 4686 




4192 4200 4208 



4534 4541 4548 



4680 4686 4692 



4325 4340 4355 
Hel 4387 




346-097 09 V 

Ha 



Hel + Hell 4026 







1.00 | 






0.80 - 






0.60 - 




1.00 






PI j 




0.80 






0.60 







6542 6562 6582 
Hel 4471 



4017 4026 4035 
Hel 4713 




4385 4388 4391 
Hell 4200 



4467 4471 4475 
Hell 4541 



4711 4713 4715 
Hell 4686 




4190 4200 4210 



4533 4541 4549 



4682 4686 4690 




346-107 09.5 V 

Ha 



Hel + Hell 4026 




4325 4340 4355 
Hel 4387 



6542 6562 6582 4017 4026 4035 

Hel 4471 Hel 4713 




t 1.00 






1.00 ; 


; 0.90 






0.95 ; 


0.80 






0.90 ; 


; 0.70 






0.85 ; 



4383 4388 4393 4463 4471 4479 4710 4713 4716 




4190 4200 4210 



4533 4541 4549 



4680 4686 4692 



Fig. A.5. Same as Fig.|AT] however, for NGC 346-090, -093, -097 and -107. 



weak Hen lines at 4200 and 4541 A could still be fitted using 
the automated method. 



NGC 346-097 The final fit to the spectrum of the 09 dwarf 
346-097 is presented in Fig. IA.5I Even though the nebular 
emission is quite severe in the hydrogen lines and in some of 
the helium lines a good fit was obtained. Its parameters are in 
agreement with other 09 V objects in our sample. Note that the 
gravity of log g = 4.5 is remarkably high. 

NGC 346-107 The nebular contamination, as shown in 
Fig. IA.5I is severe in the spectrum of this object. All hydro- 



gen and neutral helium lines suffer from nebular emission. This 
did not prevent the automated fitting method from obtaining a 
reliable fit. Most importantly, the effective temperature, deter- 
mined without information about the cores of the He i lines, is 
in good agreement with the other 09.5 V objects studied here. 



NGC 346-1 1 2 The spectrum of 346-1 12 is rather noisy, which 
is the result of its low luminosity of logL = 4.36 Lq, making 
it the intrinsically faintest object in our sample. Also strong 
nebular contamination is present in the line cores. The fit to the 
spectrum is presented in Fig. IA.6I Despite the low quality of 
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Fig. A.6. Same as Fig.ET] however, for NGC 346-1 12, NGC 330-031 and -052. 



the data the fits of the line profiles are good. The low S/N value 
does, however, result in large error bars (cf. Tab.@}. 



NGC 330-01 3 As can be seen in Fig.|A6]all lines of this 08.5 
giant could be fitted accurately. In the HR-diagram in Fig. |6] 
this object is found to lie close to the region in which rotating 
evolutionary models predict a helium surface enrichment. The 
relatively high helium abundance of Yu e = 0.18, which was 
needed to obtain this fit, seems to be in agreement with rota- 
tionally enhanced mixing effects (see Sect. l4.4> . 



NGC 330-052 Presented in Fig. IA.6I is the spectrum of 
NGC 330-052. The H5 Balmer line was not fitted as the spec- 
trum around in that region was of a too low quality. Despite 
this missing line and the low signal-to-noise ratio, a good fit 
could be obtained. The low quality of the spectrum does result 
in relatively large error estimates (cf. Tab.0J. 



AzV 1 4 For this object and the remaining objects the hydrogen 
Balmer line H/3 was also fitted. The spectrum of the 05 V star 
suffers from strong nebular contamination in both the He i and 
He n lines. However, the fit shown in Fig. IA.7I does seem to 
reproduce the observed features correctly. A comparison with 
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Fig. A.7. Same as Fig. IA.ll however, for AzV 14, 15, 26 and 95. 
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the parameters determined by Masse vet alJ J2004). who also 
used fastwind, confirm this. Within the error bars the parameter 
values are comparable. 



AzV 1 5 The He i A4-47 1 line in the spectrum shown in Fig. IA.7l 
is strongly affected by nebular emission. This hampers an accu- 
rate determination of the effective temperature and the surface 
gravity, which is reflected in relatively large error estimates 
listed in Tab.|4] Nebular emission can also introduce a system- 
atic error in the determination of these parameters. However, 
the behaviour of this object in the T e ff vs. spectral type dia- 
gram and log T e ff- log g c plane does not indicate a problem. 



Moreover, a "normal" helium abundance of Tee = 0.10 is re- 
covered from the spectrum, indicating that the He i line profiles 
contain enough information to avoid a degeneracy effect be- 
tween r e jf and Yue (also see Paper I). All in all, we conclude 
that despite the strong nebular contamination the parameters of 
AzV 15 are robust. This conclusion is strengthened by the good 
agreement between the projected rotational velocity we deter- 
mine from the optic al (135 km s~ 1 ) and the value determined by 
IPennv et"aPd2004h . These authors used cross-correlation tech- 
niques to determine v t sin i = 128 km s _1 from UV spectra. 
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Fig. A.8. Same as Fig. EOl however, for AzV 243, 372, 388 and 469. Shown with a dotted line for AzV 372 is the effect of an 
increase of the mass loss rate by 0. 1 dex on the line profiles of Ha and He n /14686. 



iHean et all §006) derive a considerably lower T e g of 
37 kK. This is connected to the log g value obtained by these 
authors that was lower by ~0.2 dex. 



AzV 26 A good fit was obtained for AzV 26, though in the fi- 
nal fit presented in Fig. lA.7l there is a small under prediction of 
the line wings of He n /14686. A possible explanation of the lat- 
ter may be an incorrect projected rotational velocity. However, 
the value u r sin i - 132 km s~ 1 obtained from the automated fit 
is in perfect agreement with the value of v, -sin / = 127 km s -1 



IPennvet al.l J2004) derived from the UV spectrum. Therefore, 
this small discrepancy remains unexplained. 



Compared to the analysis by Mas sev et al.l d2004l) the 
largest differences are found for T e g and log g, for which the 
automated method determined values higher by ~2 kK and 
0.25 dex, respectively, and M, for which a value lower by 
0.16 dex was obtained. A small part of the difference in sur- 
face gravity can be explained by the larger v T sin i of 150 km s _1 
adopted by these authors. The largest part of the difference in 
log g and we believe is the result of a better multidimen- 
sional parameter optimisation by the automated method (see 
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above). The reduced mass loss rate can be explai ned by the 
larger value of ft found in this study. Mas sev et alJ assumed a 
fixed value of/3 - 0.8, whereas we find ft = 1.17. 

AzV 95 The best fit obtained for this 07 giant is presented in 
Fig. IA.7I The accurate reproduction of all lines and the mod- 
est error estimate s require no further comments. We note that 
iHeanetal IfcOOfil) in their analysis obtained a T e ff value that is 
lower by ~3 kK. We ascribe this difference to the gravity deter- 
mination of these authors, which was lower by approximately 
0.3 dex. 

AzV 243 In Fig. lA.8l the spectrum of AzV 243 is presented. 
With the automated method we were able to obtain a very 
good fit for a ll fines. The final parameters can be compared 
to the work of Haser et al. ( 1998). With exception of the effec- 
tive t emperature that H aser et alJ adopted from the analysis of 
IPuls et all i 19961) and ft, there is good agreement. Our effective 
temperature is lower by ~1000 K, which we attribute to the in- 
clusion of line blanketing in our analysis. The acceleration of 
the wind was determined bv Hase r et alJ to be B = 0.7, based on 
UV resonance lines. Our value of ft = 1 .4 was derived from the 
optical spectrum. The reason for this discrepancy is unclear. 
However, it could be related to the fact that the wind is very 
weak (order 10~ 7 M0yr~'), introducing a large error in the de- 
termination of this parameter. 

Compared to the study o f lPennvet alJfe004 we again find 
good agreement between our v t sin i determined from the op- 
tical (59 kms -1 ) and their i> r sin/ derived from a UV analysis 
(62kms- 1 ). 



lEvans etall We can explain this difference by the large ft these 
authors derive. To match the Ha profile Evans et al.l needed 
ft = 2.25, whereas we obtained ft = 1.27 using the automated 
method. As a larger value for ft results in a slower accelerating 
wind and, consequently, in a higher density in the Ha line form- 
ing region, the required mass loss rate to reproduce the Ha line 
pr ofile is d ecreased. Note that our wind parameters, in contrast 
to Evans et all can reproduce the He n /14686 line profile cor- 
rectly. However, one should realise that in our fitting method 
th is line was g i ven a r elatively high weight compared to the fit 
of E vans et allj2004al) . who gave priority to Ha. Consequently, 
the differences in wind parameters for this object also reflect 
the particularly peculiar shape of the wind lines and the diffi- 
culties involved in reproducing them. 

AzV 388 For the 04 dwarf AzV 388 a good fit was obtained. 
As is shown in Fig. lA.8l all line profiles are reproduced accu- 
rately. The projected rotational velocity of 163 km s _I required 
to fit the profi les is in good agreement with the UV determina- 
tion by|P£nn^^0<l2QQ3)- They find v r sin i — 179kms~ 1 . 

AzV 469 The fit to the spectrum of AzV 469 is shown in 
Fig. IA.8l For all lines a good fit was obtain ed. Our fit par ame- 
ters compare well with the results of Evans et alJ (|2004a) who 
studied this object using cmfgen. In contrast we find large dif- 
ferences for log q, 7Vff and M wh en comparing our results to the 
analysis of Massev et alJ (|2Q04) who also analysed AzV 469 
using fastwind. Note, however, that the differe nces a re likely 
the result of a not well constrained fit by Massev *et all as these 
authors had difficulties matching the Ha core emission and 
severely under predicted the strength of the He i A447 1 line. 



AzV 372 Based on its peculiar Ha profile and the fact that 
they were una ble to simultaneously fit the He i and He n lines 
Mas sev et all J2004I) suggest that this object is a spectroscopic 
binary. However, the fit presented in Fig. lA.8l shows that a si- 
multaneous fit of both the He i and He n lines for a single set 
of paramet ers is possible. Therefore, we do not concur with the 
findings of iMassev et all and treat AzV 372 as a single star. 
Note that we did encounter a mild form of the "generalised di- 
lution" effect in the He i AAAll line. 

Using our fitting method, reproduction of the Ha line still 
remains difficult. We are not able to satisfactorily match the ab- 
sorption and the weak emission simultaneously. To fit the emis- 
sion an increase of M of the same order as the error estimate of 
this parameter is necessary. In Fig. lA.8l the effect of an increase 
of the mass loss rate with 0. 1 dex on the line profiles of Ha and 
He n /14686 is shown using a dotted line. Consequently, given 
the modest increase in M required we regard our determination 
of the mass loss rate as correct within the given error estimates. 

AzV 372 was also analysed by Evans etafl J2004a). 
Compared to this study we find differences for T e g, Yu e , M 
and ft. The effective temperature from the automated method 
is larger by approximately 3 kK. This is probably linked to the 
reduced helium abundance we find and t he improved fit of the 
He n lines (e.g. see Fig. 3 of IE vans et all) . Our mass loss rate is 
approximately a factor of two larger than the M determined by 
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